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METHOD FOR NONINVASIVE CONTINUOUS 
ETERMINATION OF PHYSIOLOGIC CHARACTERISTICS 

lELD OF THE INVENTION 

The present invention relates generally to noninvasive methods of quantitatively 
determining various physiologic parameters relating to cardiovascular and respiratory 
function. More particularly, the invention relates to methods and devices for continuous, 
noninvasive determination of: hemoglobin, glucose and other blood constituent 
concentrations, blood pH and acid-base balance, blood flow differentials, blood 
temperature, arterial blood pressure, venous pressure, arterial oxygen saturation, venous 
oxygen saturation, arterial pulse wave velocity, aortic pulse wave velocity, aortic pulse 
flow velocity, cardiac stroke volume, cardiac index, cardiac output, heart rate, respiratory 
rate and cardiac ejection fraction. 
BACKGROUND OF THE INVENTION 

Critically ill and seriously injured patients require constant care and attention. 
Doctors, nurses, and hospital technicians need a continuous flow of information about the 
many patients under their care. Heart rate and blood pressure measurements are two 
primary vital signs that indicate the health of patients under their care. When these two 
common indices of wellness fall below normal readings, a patient is usually in distress and 
requires immediate attention. 

Dangerous conditions brought about by a cardio-vascular or pulmonary disease, 
severe trauma, or drug abuse may bring about a failure of the lungs and heart to supply the 
bloodstream with life-giving oxygen. Such a fatal deficiency can be detected by 
continually gauging the amount of hemoglobin in the bloodstream that is carrying oxygen. 
This third vital sign, which manifests oxygen saturation of the blood, is especially critical 
because a rapid decline in oxygen in the bloodstream is associated with increased risk of 
patient mortality. 

It is well known that blood pressure can be directly measured by placing a fluid- 
filled catheter directly into the vessel and coupling this to an electro-mechanical 
transducer. This is the most accurate means, but has all the disadvantages of invasive 
measurement, including pain on insertion, risk of infection or disease transmission, risk of 
bleeding or thrombosis, and great expense. A further disadvantage is the creation of toxic 
medical waste (needle, gloves, skin dressing, etc). 
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Blood pressure measurement can also be measured indirectly using an occlusive 
cuff (with either auscultation or oscillometry to make the determination). This is the most 
common means of blood pressure measurement. Illustrative is Pat. Nos. 5,582,179, 
5,048,533, 5,152,296 and 4,793,360. 

A further occlusive cuff apparatus is disclosed in U.S. Pat. No. 5,766,130. 
According to the invention, the apparatus includes multiple "pressurized pneumatic cuffs" 
that are used to "plot blood pressure and/or volumetric blood flow wave forms from a 
plurality of separate digits and/or extremities of a patient so that circulatory parameters 
may be measured rapidly and recorded from a great number of the patient's digits or 



Although commonly employed, the occlusive cuff also has numerous 
disadvantages, which include discomfort, intermittent readings, and poor reliability. 

An additional means of determining blood pressure is through an assessment of 
"pulse wave velocity". Several prior art references disclose methods and/or apparatus 
employing such means. Illustrative is U.S. Pat. No. 5,649,543. 

There are also several prior art references that disclose methods and/or apparatus 
for determining blood pressure through a "pulse wave amplitude" assessment. Illustrative 
are U.S. Pat. Nos. 4,735,213, 4,872,461, 4,793,360, 5,265,01 1, 5,385,149, 5,511,303, 
5,582,179, 5,680,867 and 5,882,31 1. 

Additional physiologic characteristics such as blood temperature and pH provide 
further information regeirding the status of the patient. Moreover, combinations of 
measurements can be used to determine specific cardio-pulmonary parameters. 

Acid-base balance (the most common measure is pH) is perhaps the most 
important factor in the chemistry of both biologic and non-biologic systems. It figures in 
speed of reactions; indeed, if a reaction will occur at all. In most biologic systems, 
determination of pH requires laboratory analysis. The monetary costs are high, and 
procedures involve risk for patient subject and laboratory technicians among others. Toxic 
medical waste (syringes, gloves, etc.) is created and must be disposed of safely. In other 
systems (and to some extent in biologic systems), pH measurement is done in one of two 
common ways: colorimetric and electrochemical. Noninvasive measurement of arterial 
blood pH is described in U.S. Patent No. 5,978,691. 

Although most of the noted noninvasive monitoring methods and apparatus, 
particularly the occlusive cuff, have been employed for many years by health care 
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personnel, the conventional methods and apparatus have one major, common drawback - 
the need for separate calibration. 

Accordingly, there is a need for noninvasive methods and devices capable of 
continuously determining various physiological characteristics, such as blood pressure, 
central venous pressure and cardiac output, without separate calibration. There is also a 
similar need for noninvasive methods and devices for determining various blood 
parameters including hemoglobin, glucose and other blood constituent concentrations, 
blood pH and acid-base balance, blood flow differentials, blood temperature, blood 
pressures and pressure wave differentials. As will be appreciated by one having ordinary 
skill in the art, the present invention satisfies these and other needs. 
SUMMARY OF THE INVENTION 

The invention comprises methods for noninvasively determining the concentration 
of a blood constituent, generally including the steps of providing a tissue probe having a 
radiation emitter with a wavelength and a detector configured to receive the wavelength 
after absorbance through a path length of the patient's blood; measuring absorbance of the 
patient's blood by emitting and detecting radiation after passage through the patient's 
blood; varying the path length of to provide multiples of path length; measuring 
absorbance of the patient's blood at each multiple of the path length; and determining the 
concentration of the blood constituent based upon the changing absorbance. Preferably, 
the determination is made continuously. The determination can also be made using 
calculated multiples of the path length. 

Alternatively, the invention comprises the steps of providing a first and second 
tissue probe each having a radiation emitter with a wavelength and a radiation detector 
configured to receive the wavelength after absorbance through a path length of the 
patient's blood at a position relative to the heart of the patient; measuring absorbance of 
the patient's blood by emitting radiation and detecting radiation after passage through a 
first path length of the patient's blood; varying the pressure of the blood within the first 
and second probes; measuring absorbance of the blood as the pressure is varied; and 
computing the time of arrival and amplitude of the pulse based on the absorbance at the 
varying pressures. Preferably, the pressure of the blood within the probes can be varied by 
changing the hydrostatic pressure relative to the heart such as by changing the height of 
the probe. 
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These techniques can be applied to arterial or venous blood. Preferably, the blood 
constituent being measured comprises hemoglobin. By measuring hemoglobin oxygen 
saturation or blood pH at different probe locations, blood corresponding to a single flow 
wave can be identified. By comparing the timing of the flow wave arrival at the different 
probe locations, multiple determinations of cardiac characteristics can be made, including 
arterial pulse wave velocity, aortic pulse wave velocity, aortic pulse flow velocity, cardiac 
stroke volume, cardiac index, cardiac output, heart rate, respiratory rate and cardiac 
ejection fraction. In other embodiments, tissue probes having more than one radiation 
emitter and detector pairs can be employed. 

The invention also comprises methods for noninvasively determining the pH of 
blood of a patient, generally including the steps of providing a first tissue probe having a 
first and second radiation emitter with a first and second wavelength, respectively, and a 
first and second radiation detector configured to receive the first and second wavelength, 
respectively, after absorbance through the patient's blood; measuring absorbance of the 
patient's blood by emitting radiation at the first and second wavelength through the 
patient's blood and detecting the radiation after passage through the patient's blood; and 
computing the pH of the blood based upon the measured absorbance at the first and second 
wavelengths. The first and second wavelengths are selected so that the absorbance of the 
first wavelength depends upon the pH of the blood and the absorbance of the second 
wavelength is substantially independent of the pH of the blood. Preferably, the method is 
tailored to measuring the absorbance spectrum of hemoglobin species that have certain pH 
dependent absorption peaks. Since pH is also dependant upon temperature, more accurate 
results can be obtained if the blood temperature is varied. Alternatively, if the pH of the 
blood is characterized, the temperature of the blood can accurately be determined. 

Yet another embodiment of the invention is a method for non-invasively 
determining the concentration of a blood constituent comprising the steps of measuring 
absorbance of arterial and venous blood; determining arterial and venous oxygen 
saturation; subtracting hemoglobin absorbance based upon the arterial and venous 
saturation; and determining the concentration of a blood constituent based upon remaining 
absorbance. Preferably, the blood constituent^comprises glucose. 

The invention also includes a device for the noninvasive monitoring of a 
physiologic characteristic of a patient's blood. In one embodiment, the device comprises a 
tissue probe having a radiation emitter and a radiation detector configured to receive the 
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radiation after absorbance through the patient's blood; a position sensor for determining 
the relative height of the probe compared to a level corresponding to the patient's heart; 
and a controller for computing the physiologic characteristic of the patient's blood based 
on the absorbance of the first wavelength of radiation and the relative height of the probe. 
The radiation emitters of the invention can utilize a single wavelength or a plurality of 
discrete wavelengths and may include visible light, infrared light, and ultraviolet light. 
The probes are adapted for use with hands, fingers, feet, toes, ears, earlobes, nares, lips, 
tongue and the like. Additional radiation emitters and detectors may also be used. 
Preferably, the probe further comprises ECG leads. 

An alternative embodiment of the device of the invention comprises a tissue probe 
and controller in conjunction with a movement generator for inducing a position change of 
the probe with respect to a level corresponding to the patient's heart. Preferably, the 
movement generator induces a known position change of the probe and moves the probe 
to positions above and below a level corresponding to the patient's heart. 

The invention also comprises methods for determining other physiological 
characteristics of a patient's blood noninvasively. In one embodiment, absorbance 
characteristics of the blood are measured at varying positions relatively to the level of the 
patient's heart. By comparing blood parameters such as pulse amplitude, pulse velocity, 
pulse delay, pulse contour, flow velocity and flow delay to hydrostatic pressure 
differences induced by the position changes, characteristics such as arterial and central 
venous blood pressure and cardiac output can be determined. Alternatively, two probes 
are used to compute pulse delays between coupled tissues or opposing tissues. 

The subject invention relates novel methods for noninvasive determination of 
physiologic characteristics. The first new and unique method and device utilizes changes 
in hydrostatic pressure induced by positional changes to facilitate measurements. A 
second new and unique method and device for noninvasive determination of cardiac 
output by measuring delays in pulse arrival times in coupled organs or members on 
opposite sides of the body is also described. The two methods are such that they can 
advantageously be used together. 

By varying the hydrostatic pressure in an extremity, one can not only perform self- 
calibration for a blood pressure determination, but also change the pulse wave velocity and 
pulse propagation delay with respect to the opposite extremity. With this information. 
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pulse wave velocity, and consequently flow wave velocity at the aortic root can be 
determined. 

Similar techniques of varying hydrostatic pressure can be used to assess venous 
pressure and saturation. The technique of repetitious determinations made while altering 
position or other variables allows a multitude of additional analyses to be made. The 
determinations can be made intermittently or continuously. 

Further objects of the invention are exemplified by the following potential 
applications: 

(A-1). A patient is anesthetized for a surgical procedure. Probes are attached to 
the index fingers of each hand, and a movement generator is placed on one arm. A 
complete set of vital signs and physiologic characteristics is generated continuously, 
including: arterial blood pressure, venous pressure, arterial oxygen saturation, venous 
oxygen saturation, arterial pulse wave velocity, aortic pulse wave velocity, aortic pulse 
flow velocity, cardiac stroke volume, cardiac output, heart rate, and respiratory rate. Other 
chziracteristics can be calculated if desired. 

(A-2). A patient is anesthetized for a cardiac surgical procedure. As access to the 
arms is difficult, probes are attached to the patient's temples. A complete set of vital signs 
and physiologic characteristics is continuously generated. 

(A-3). A patient is anesthetized for a cardiac surgical procedure; this time the 
procedure includes valvular repair or replacement. Since the cardiac output and other 
characteristics can be continuously computed, the adequacy of the surgical repair can be 
judged immediately. 

(A-4), As the number of endoscopic or minimedly invasive cardiac surgical 
procedures is expected to increase, the demand for less invasive monitoring will also 
increase. The device described herein provides noninvasive, continuous monitoring of 
essentially all cardiovascular characteristics. 

(A-5). Cardiac catheterization procedures are often done on critically ill patients. 
As the procedures are usually relatively brief and accomplished without general 
anesthesia, invasive monitoring methods are often not desired despite the illness of the 
patients. The device described herein will provide the necessary monitoring that is 
typically provided by much more invasive, expensive, and time consuming monitors 

(A-6). A patient is hospitalized in the intensive care unit of a hospital after a heart 
attack. Probes are attached to the index fingers of each hand, and a movement generator is 
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placed on an arm or a leg. A complete set of vital signs and physiologic chEiracteristics 
can be continuously generated. In addition, arrhythmias can be detected and diagnosed. 

(A-7). The patient noted above is now moved to a "step-down" or telemetry unit 
from the intensive care unit. Because the device described herein eliminates the need for 
5 invasive monitoring lines, a complete set of vital signs and physiologic characteristics can 
still be continuously generated. As the patient has mobility of arms and legs, a movement 
generator is no longer needed, as the patient's spontaneous motion, even during sleep, will 
generate hydrostatic pressures in the limbs, allowing all computations to be made. In 
addition, the probes may be made wireless, and connected to a central nursing station by 
10 means of infrared or radio frequency communication. 

(A-8). The patient noted in applications 6 and 7 above is now moved to a regular 
hospital bed, and does not require continuous monitoring. However, vital signs can still be 
recorded by a technician moving the device from bedside to bedside on a cart. The device 
P does not require highly trained nursing personnel to operate. 

Q 1 5 (A-9). The patient noted in applications 6, 7, and 8 above has now been 

?S discharged from the hospital, and now presents to his physician's office for follow-up. 

The same device can be used in physician's offices, as it provides better care at lower cost. 
(A- 10). Ambulances, emergency vehicles, and military vehicles can also employ 
1=* this device as it is very simple to operate, and provides data that currently is impossible for 

20 them to obtain. In addition, the information can be transmitted to central stations where 
l'^ medical personnel are available for help and advice. 

|:i (A-1 1). The device and methods of the invention will provide means of 

monitoring patients or checking vital signs for extended care facilities, nursing homes, and 
other health-related facilities 

25 (A- 12). Blood pressure screening clinics and drugstores will have a greatly 

improved means of determining patient's blood pressures and other vital signs. Airports 
and airplanes are able to purchase medical equipment, but often do not have personnel 
trained to operate the equipment. The device is simple and quick to operate. 

(A- 13). The patient noted in applications 6 through 9 above can also monitor his 

30 heart disease and health care at home. The operation of the device is straightforward 

enough to be used by the layman with minimal instruction, and inexpensive enough for 
personal home use. The patient can measure his cardiovascular characteristics daily, or as 
frequently as he and his physician desire. A communication means, such as a modem, can 
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easily be incorporated into the device. This, with appropriate software and support, would 
allow essentially instantaneous communication with a physician's office, clinic, or 
hospital. In addition, a permanent record can be made and stored electronically. If 
desired, the device could automatically "sign on" to the Internet or other network, and link 
to the appropriate website or other address. The ability to participate more fully in their 
own health care will improve the welfare of individuals. 

(A-1 4) The patient of above presents to the emergency room of a hospital wdth 
chest pain. The ER physician can access, via the Internet or other means, the patient's 
vital sign history, including ECG. This allows the physician to determine if abnormalities 
are new or chronic. Changes, such as dysrhythmieis, can be identified as to when they first 
occurred, perhaps to within a time frame of hours or less. 

(A-1 5). People without diagnosed cardiovascular disease can use the device to 
allow themselves to participate in their own health care. This will allow virtually 
immediate diagnosis of any problems, allowing early intervention. In addition, a 
permanent record can be created if desired. 

(A- 16). The device will impact fitness and physical training for everyone from lay 
people to military personnel to professional athletes. 

(A- 17). The device can be employed in the diagnosis and management of 
peripheral vascular disease. Measurement of pulse wave velocity in the extremities, and 
particular differential pulse wave velocities in the lower extremities, can be used to 
diagnose peripheral vascular disezise. Since measurements are real time and continuous, 
they can also be used in management. For example, if balloon angioplasty of an artery is 
performed, the clinician can tell immediately if flow has improved. In the case of 
angioplasty of coronary arteries, the clinician can follow cardiac characteristics on a beat- 
by-beat basis. 

(A-1 8). In addition to peripheral vascular disease, other diseases, such as 
abdominal aortic aneurysm, can be diagnosed and managed. Changes in pulse wave 
velocity and waveform can be followed for years if desired. 

(A-1 9). Some of the most important potential uses of the device relate to the 
health care of neonates and young children. For these patients, the measurement of 
common characteristics such as blood pressure can be difficult even for highly trained 
personnel in well-equipped facilities. The simple placement of probes on fingers will 
alleviate this. The device will also allow noninvasive diagnosis of congenital cardiac 
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defects and anomalies. Analysis of differential pulse wave velocity and blood pressure 
will allow rapid, accurate, and specific diagnosis of many disorders, including Tetralogy 
of Fallot and transposition of the great vessels. The ability to distinguish both arterial and 
venous saturations and pressures will allow diagnosis of patent ductus arteriosus, truncus 
arteriosus, atrial septal defect, and ventricular septal defect. Differential arm and leg pulse 
wave velocities and pressures will confirm diagnosis of coarctation of the aorta. Because 
of its continuous measurements, the device can be used for only for diagnosis but 
confirmation of adequacy of repair, including intraoperatively. As the device is 
inexpensive and easy to operate, it may become a screening tool for newborns and infants. 

(A-20). The device can be used in conjunction with intra-aortic balloon pump 
(lABP) counterpulsation. Beat-by-beat analysis of effectiveness and ability to wean fi-om 
counterpulsation can be made. 

(A-21). The device can be used in conjunction with placement of cardiac 
pacemakers, to set proper rate and timing intervals. In addition, efficacy of pacemakers 
can be checked as frequently as desired, and scheduling of reprogramming or replacement 
made automatically. 

(A-22). It is straightforward to incorporate other devices, such as the 
electroencephalogram (EEG) or electromyogram (EMG), into probes of the invention. As 
a general-purpose monitor, the device will invite the addition of specialized add-ons. 

(A-23). Many enhancements are included in the invention. For example, addition 
of chest (horizontal) leads allows full diagnostic ECGs to be performed. 

(A-24). Under some circumstances, such as severe hypotension, the pulse cannot 
be identified in the periphery. In such cases, many of the determinations claimed herein 
cannot be made. However, the ability of the device to identify venous blood can still give 
important information. 

(A-25). Forces other than gravity can be used. In a microgravity environment 
such as a space station orbiting the Earth, a device such as the one described could be 
constructed to perform all indicated determinations using acceleration caused by 
movement in place of gravitational acceleration. 

(A-26). As mentioned in the examples above, an anticipated use is in the field of 
home health care, with the possibility of automatic sign-on and direction to a website. As 
the user is already participating in his or her health care, the extension of providing access 
to related health or other information via the Internet® is a natural one. 
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(A-27). A verification means, such as fingerprint scanning, can be incorporated 
into a personal-use device, to ensure that any medical information gathered belonged to 
the individual using the device. 

(A-28). The device will be used in conjunction with the Penaz technique or other 
methods, such as calibration with a cuff or other means, as desired. 

(A-29). Nerve blocks, such as spinal or epidural anesthesia, block the actions of 
the sympathetic nervous system, producing both arterial and venous dilation. This device 
can be used to assess these effects. 

(A-30). A patient with suspected peripheral vascular disease of the lower 
extremities is tested. A probe of the invention is placed on a toe of one foot, and an 
identical probe is placed on a toe of the other foot. Pulse velocities and possible pulse 
delay between the two extremities aire meeisured. 

(A-3 1). The patient fi*om application 30 is tested further. Another identical 
probe is placed on a finger. Pulse velocities and possible pulse delay between the two 
lower extremities are measured. In addition, comparison is made with pulse velocity and 
propagation time in the upper extremity. 

(A-32). The patient from applications 30 and 31 is tested further. This time, the 
legs are raised and lowered to gain further information on pulse velocities and propagation 
times. 

(A-33). A person suffers severe traimia to the lower extremities in a motor vehicle 
accident. His physicians are concerned about the possibility of his developing 
"compartment syndrome", which is compromise of blood supply due to swelling of 
injured tissues. Probes of the invention are placed on his feet to monitor the integrity of 
the vascular supply. 

(A-34). A pharmaceutical company develops a new drug to manage hypertension, 
CHF, arrhythmias. Instead of long trials or invasive trials, they are conducted on subjects 
using this device. 

(A-3 5). Taking vital signs on a space station. The techniques described herein 
lend themselves well to the special problems of this unique environment. 
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BRIEF DESCRIPTION OF THE FIGURES 

Further features and advantages will become apparent from the following and more 
particular description of the preferred embodiments of the invention, as illustrated in the 
accompanying drawings, and in which like referenced characters generally refer to the 
same parts or elements throughout the views, and in which: 

FIGURE 1 is a diagram of the central cardiovascular system, showing the 
asymmetry of origins of the vessels off the aortic arch; 

FIGURE 2 is a partial perspective view of one embodiment of a physiologic 
"tissue" probe with a single emitter-detector pair, according to the invention; 

FIGURE 3 is a partial perspective view of an alternative embodiment of a tissue 
probe with a single emitter-detector pair, according to the invention; 

FIGURE 4 is a partial perspective view of an alternative embodiment of a tissue 
probe with two emitter-detectors, according to the invention; 

FIGURE 5 is a partial perspective view of an alternative embodiment of a tissue 
probe with a single emitter-detector pair and a single electrocardiogram (ECG) electrode, 
according to the invention; 

FIGURE 6 is a partial perspective view of an alternative embodiment of a tissue 
probe with a single emitter-detector pair and two ECG electrodes, according to the 
invention; 

FIGURE 7 is a partial perspective view of an alternative embodiment of a tissue 
probe with a two emitter-detector pairs and two ECG electrodes, according to the 
invention; 

FIGURE 8 is a partial perspective view of an alternative embodiment of a tissue 
probe having a position sensor, according to the invention; 

FIGURE 9 is a schematic illustration of a patient with tissue probes placed on 
opposite digits according to the invention; 

FIGURE 1 0 is a schematic illustration of a patient with probes placed on opposite 
temples, according to the invention; 

FIGURE 1 1 is one embodiment of a circuit schematic including a 
photoplethysmogram, according to the invention; 

FIGURE 12 is an alternative embodiment of a circuit schematic having a 
photoplethysmogram with an ECG amplifier, according to the invention; 
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FIGURE 1 3 is an alternative embodiment of a circuit schematic having a 
photoplethysmogram with an ECG amplifier and a level signal, according to the invention; 

FIGURE 14 is an alternative circuit schematic having a photoplethysmogram with 
two independent channels, according to the invention; 

FIGURE 1 5 is an alternative circuit schematic having a photoplethysmogram with 
two independent channels and an ECG amplifier, according to the invention; 

FIGURE 1 6 is a schematic illustration of a patient with tissue probes placed on the 
digit and arm near the brachial artery, according to the invention; 

FIGURE 17 is a schematic illustration of a patient with tissue probes placed on a 
finger and toe, according to the invention; 

FIGURE 1 8 is a schematic illustration of a patient with tissue probes placed on 
opposite fmgers and a toe, according to the invention; 

FIGURES 1 9 and 20 are schematic illustrations of patients with tissue probes 
placed on opposite digits and positioned at differential heights relative to the patients' 
heart, according to the invention; 

FIGURE 21 is a schematic illustration of a patient with tissue probes placed on a 
toe and opposite fingers positioned at different heights, according to the invention; 

FIGURES 22-25 are graphical illustrations of oscilloscope tracings showing 
recordings using methods of the invention wherein the probes are placed on index fingers 
of opposite hemds, and the pulse velocities and relative arrival times are varied by 
changing the heights of the probes relative to the heart; 

FIGURE 26 is a graphical illustration of oscilloscope tracings showing trains of 
pulses with varying oxygen saturation values; 

FIGURE 27 is a graphical illustration of oscilioscope tracings showing trains of 
pulses with varying pH values; 

FIGURE 28 is a graphical illustration of oxygen saturation versus heart rate 
showing the effects of induced changes in oxygen saturation measured at the toe; 

FIGURE 29 is a graphical illustration of oxygen saturation versus heart rate 
showing the effects of induced changes in oxygen saturation measured at the left hand; 

FIGURE 30 is graphical illustration of oxygen saturation versus heart rate showing 
the effects of induced changes in oxygen saturation measured at the right hand; 
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FIGURES 31 and 32 are graphical illustrations of oxygen saturation (SO2) versus 
pH showing the effect of plasma pH on oxygen saturation mezisurements for umbilical 
cord blood samples equilibrated with 100% oxygen; 

FIGURE 33 is a computer algorithm for determining hemoglobin concentration, 
according to the invention; 

FIGURE 34 is a graphical illustration of optical extinction curves of four species 
of hemoglobin at light wavelengths in the range of 600 to 950 nm; 

FIGURE 35 is a graphical illustration of optical extinction curves for hemoglobin 
A at four different pH values and at light wavelengths in the range of 500 to 600 nm; 

FIGURE 36 is an absorption spectrum for oxyhemoglobin between measurements 
at pH 7.6 and pH 6.8; 

FIGURE 37 is an absorption spectrum for carboxyhemoglobin between 
measurements at pH 7.6 and pH 6.8; 

FIGURE 38 is a graphical illustration of a normal hemoglobin-oxygen dissociation 
curve (HODC); 

FIGURE 39 is a graphical illustration of HODC s showing shifts due to pH and 
temperature; 

FIGURE 40 is a graphical illustration of oxygen saturation versus heart rate 
showing the change of oxygen saturation corresponding to high cardiac fraction ejection; 
and 

FIGURE 41 is a graphical illustration of oxygen saturation versus heart rate 
showing a change of oxygen saturation corresponding to low cardiac fraction ejection. 
DETAILED DESCRIPTION OF THE INVENTION 

Functionally the heart is divided into two sides or sections. The right or pulmonary 

circulation section that receives blood from the veins of the body and pumps it through the 
lungs and the left or systemic circulation section that receives the blood from the lungs and 
pumps it to the body. The blood is then collected in the veins to be returned to the right 
side of the heart. 

Referring to Figure 1 , the arterial system begins at the aorta 1 , to which the left 
ventricle of the heart pumps. The first three branches of the aorta are the brachiocephalic 
or innominate artery 2, the left (common) carotid artery 3, and the left subclavian artery 4. 
The brachiocephalic artery branches into the right subclavian 5 and right (common) 6 
carotid arteries. These arteries provide the blood supply for the head and upper 
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extremities. The aorta then passes down (caudal) through the body, continuing to provide 
arterial branches to organs, terminating as a bifurcation creating the iliac arteries. The 
brachiocephalic or innominate artery is the first branch of the aorta. It in turn branches 
into the right subclavian and right carotid arteries. In contrast, the left subclavian and left 
carotid arteries originate directly off the aortic arch. Thus, the subclavian and carotid 
arteries and £iny of their branches to one side of the body have different paths from their 
counterparts on the opposite side of the body. 

Because of the different origins from the aorta and different branching pattern of 
the arterial tree, it can be appreciated that blood ejected from the left ventricle will not 
follow symmetrical pathways to opposite arms or body or opposite sides of the head. 
Similarly, the pressure pulse wave associated with left ventricular ejection will follow 
different pathways, and can be expected to arrive at different times for paired organs or 
members of the upper body, or members of the upper body compared to the lower body. 

Measurements performed by the inventor have shown this delay can range from 
less than one millisecond to several milliseconds, depending on the subject and 
circumstances. In addition, the inventor has found that this delay can be altered by several 
methods disclosed herein. This propagation delay, its alterations, and other factors make 
possible several determinations heretofore not possible by noninvasive means. 

Blood pressure is the pressure exerted by the blood within a vessel upon the wall of 
the vessel. It is measured in units of force per unit area. Central venous pressure is the 
pressure within the large veins in the chest and the right atrium, which is the common 
emptying point for the venous system. Cardiac output is the amount of blood pumped by 
the heart, expressed in units of volume per time. 

Central venous pressure (CVP) is defined as the distending pressure present in the 
veins in the chest (proximate to the heart), and is considered equal to the pressure in the 
right atrium (which is the emptying point for the venous system). Pressure should be the 
same throughout the venous system, but there are valves to ensure that the blood does flow 
back toward the heart (for example, when standing the venous blood must flow uphill, and 
there is no pump as on the arterial side). 

As discussed in detail below, the present invention generally includes a radiation 
emitter having at least one wavelength being applied through a patient's tissue to the 
patient's blood; a radiation detector which detects reception of the at least one wavelength 
after absorbance through the blood, a movement generator for inducing position changes 
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in the tissue; and a controller for computing the various characteristics based on the 
absorbance of the at least one wavelength of radiation at various position levels. In a 
preferred embodiment, the radiation emitter and detector are inserted in a probe which can 
be placed about the tissue/blood to be measured. A number of suitable configurations for 
probes are shown in Figures 2-8 and disclosed in Co-Pending Application S/N 09/684,104, 
which is incorporated by reference herein. 

For example, Figure 2 shows a representative probe 1 0 with a single emitter- 
detector pair 12. The emitter and detector are placed such that transmittance through a 
body member, such as a fmger 13, is measured. Generally, any part of the body that can 
be successfully transilluminated with the radiemt energy used can be utilized. Thus, toes, 
ears, etc. could also be used. In addition, pulse oximetry can be accomplished with this 
and all of the following embodiments. Figure 3 shows a representative probe 14 with a 
single emitter-detector pair 16 placed such that reflectance of a body member, such as a 
finger, is measured. Further, Figure 4 shows a probe 1 8 with two emitter-detector pairs 20 
and 22 spaced a known distance apart. This can be used to measure the velocity of the 
pulse wave within the probe itself. 

In certain embodiments of the invention, the probe comprises one or more 
electrocardiogram (ECG) electrodes in conjunction with the emitter-detector pairs. For 
example. Figure 5 shows a probe 24 with a single emitter-detector pair 26 and a single 
electrocardiogram (ECG) electrode 28. Similarly, Figure 6 shows a probe 30 with a single 
emitter-detector pair 32 and two ECG electrodes 34 and 36 and Figure 7 shows a probe 38 
with two emitter-detector pairs 40 and 42 and two ECG electrodes 44 and 46. Such 
probes, if placed on opposite extremities of a patient, can be used to measure central and 
peripheral pulse wave velocity as well as ECG. Other configurations, such as double 
emitter-detector pairs and single ECG electrode, can be envisioned. 

In yet other embodiments of the invention, the probe further comprises a position 
sensing or measuring device together with the emitter-detector pairs and/or ECG 
electrodes. Figure 8 shows a probe 48 similar to that shown in Figure 7 with the addition 
of a position sensor 50. This position sensor could be used in conjunction with a position 
sensor placed at heart level in order to determine the hydrostatic pressure difference 
between the two position sensors. 

As discussed in detail herein, the invention employs hydrostatic pressure to enable 
precise self-calibration of the devices in a completely noninvasive manner. Hydrostatic 
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pressure affects all liquids. Gravity or other acceleration will affect both the arterial and 
venous sides of the circulation. It affects all aspects of the blood pressure equally - mean, 
systolic, diastolic. For example, an increase in height that causes a change of 1 0 torr will 
change every pressure measurement during the cardiac cycle by this amount. 
5 For example, if the "true" blood pressure (taken level with the heart) is 120/80, 

when the arm is raised an amoimt needed to decrease the measured pressure by 1 0 torr, the 
measured pressure in the arm will be 11 0/70. The pulse pressure will be the same, but the 
transmural pressure will be 10 torr lower at all times. In addition, the vessel will be 
smaller at all points. 

10 The heart is taken to be the center of the circulatory system, and all values are in 

reference to it. This is not necessary for the practice of the invention, but serves as 
reference points for values in the current medical literature. 

The electromagnetic radiation in this description will refer to light in the visible 

C3 and infi-ared range although, as noted in the attached claims, it is conceivable that other 

£5 

1 5 forms could be used. 

^2 Similarly, while the present invention primarily describes the use of 

transillumination, it will be appreciated that reflectance spectrophotometry may 

in 

alternatively be employed. 
H Operating Principles 

p 20 It is well known that Incident radiation passing through a body part is attenuated 

„^ (absorbed) in the tissue. The theoretical basis for spectrophotometic techniques is Beer's 

Q 

i=« law (the Beer-Lambert-Bouguer law) which expresses the incident intensity in terms of 

transmitted intensity and extinction coefficients of the tissue compartments through which 
the radiation has passed. The equation can be written as: 



25 Eq. 1 ln(l / lo) = E*C*L 

where: 

lo = the incident intensity of the source radiation; 

I = the transmitted intensity of the source through the sample; 

E = the extinction coefficient of the component of interest; 

30 C = the concentration of the component in the tissue itself; 

L = the optical path length (distance) through the absorber; and 



E*C*L = absorbEmce, 
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Beer's law and the practice of spectrophotometry and oximetry have been 
exhaustively reviewed in the literature. Generally, pulse oximetry in effect filters out 
signals other that pulsating (AC). In the body, it can be assumed that the pulsatile 
component of the signal is arterial blood, while all other tissue absorbers should be non- 
5 pulsatile (DC). 

An additional feature of this invention, not found in any previous disclosure, is the 
use of hydrostatic pressure changes to vary not only the amount of arterial blood, but also 
the amount of venous blood within a body member such as a fmger. Thus, hydrostatic 
changes can be used in a similar manner to the pulse to perform measurements on both 
10 arterial and venous blood. If a finger is contained within a probe, raising the probe will 
lower the hydrostatic pressure of all vessels in the finger, both arterial and venous. Both 
arteries and veins (and arterioles and venules) vydll be smaller due to lower pressure 
distending their walls. Most change will occur on the venous side of the circulation due to 

□ lower pressure. Total absorbance of the finger will decrease. As the arterial oxygen 

P 

p 15 saturation can be measured by pulse oximetry, the venous oxygen saturation can be 

W calculated in a similar manner. 

"J 

A light signal of a known intensity and wavelength is produced by means of light- 

in 

emitting diodes (LEDs) as in currently used oximeters or, as in one possible embodiment, 
1*" a broad-band light source whereby wavelengths are isolated by a rotating filter or diffusion 

□ 20 grating. In the latter case, the emitted light is distilled through a filter which allows a 

known wavelength and intensity of light to penetrate. Use of tunable lasers or other 
M equipment is also possible. If the light source is proximate to the point of use, no further 

mode of transmission will be needed. If it is not, the light will be transported to the 
desired point by mezms such as a fiber optic cable, preserving the wavelength and 
25 intensity. 

Several means of motion induction are possible. While a movement induction 
means is described herein, voluntary movement or other means are certainly possible. 
Various means of position measurement are also possible. For example, a liquid filled 
tube with an end open to the atmosphere can be employed. Other position sensors are 
30 known to those having skill in the art, and include electromagnetic, spectroscopic, 

ultrasonic, and chemical means. Preferably, a broad-band photo detector (in the case of 
visible or infrared light) or other means is employed to measure the quantity of transmitted 
light. 
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To generate a single data point, the movement induction means is used to bring the 
finger (or other space of interest) to a known position relative to the heart. Light of known 
wavelength and intensity is emitted (and transmitted if necessary) on the surface of 
interest. Detection of the light signal at a distinct point (normally opposing surface) is 
made and the relative absorbance and extinction of the signal is calculated. Signal 
processing is used to determine the pulsatile portion of the signal. The arrival time of the 
pulse is recorded, as is the amplitude and waveform. This measurement may be repeated 
one or more times to ensure the accuracy of the measurement; this can be done within a 
very short time frame (less than a millisecond). 

To generate multiple data points, the process outlined in the previous step vsdll be 
repeated at the next chosen wavelength, while still at the same predetermined position. 
The range and number of wavelengths can be selected, and changed for different 
applications. 

Once the desired number of wavelengths has been examined, the movement 
induction means would bring the finger or other volume to a predetermined second 
position, and the data collection of steps would be repeated. At the completion of 
measurements and determinations for this second position, the movement induction means 
will bring the space to a third predetermined position, and the measurements and 
determinations repeated. This process would be continued until the desired range of 
positions has been scrutinized. 

In order to make computations of pulse propagation delay, identical measurements 
would be made simultaneously with a probe on the same member on the opposite side of 
the body. For example, if one probe were placed on the index finger of the right hand, the 
other probe would be placed on the index finger of the left hand. 

Because the arterial path to the arm is essentially identical after the second part of 
the subclavian artery, any differences in pulse wave velocity and pulse wave propagation 
time must occur prior to this point; that is, very close to the root of the aorta. In any case, 
pulse wave velocity increases rapidly as the pulse wave propagates down the aorta and 
into the periphery (Fung). Thus, any timing differences in the periphery will be greatly 
reduced by the high wave velocity, leaving central effects as the most prominent. 

The apparatus of the invention can be operated intermittently or continuously. In 
the intermittent mode, a single set of calculations can be used for analysis to produce the 
determinations claimed. However, the device can also be easily operated in continuous 
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mode, with the process outlined above repeated as often as wished (constantly if desired). 
In addition, a rapid ("stat") mode can be offered with the minimum number of 
measurements made that will provide an accurate estimation of correct values. Such a 
rapid mode would be useful in emergency situations. 

While this methodology should give precise values, further adjustment may be 
desired to compensate for any discrepancies between theoretical and in vivo 
measurements. Contemporary oximeters in fact use a calibration curve when determining 
oxygen saturation, with the curve being generated with data from normal volunteers. 
Calculations and Analysis 

The following algorithms are further examples of the use of the present invention. 
Some variables have degrees of co-dependence. In these cases, values are calculated by 
iterative computational techniques. 

According to the invention, measurement of pulse wave amplitude and timing is 
made using probes such as that shown in Figure 2, using methods similar to standard 
oximetry described in the prior art. Referring to Figure 9, a first probe 52 is placed on a 
fmger and set at a known position relative to the heart. Another, simultaneous 
measurement of pulse wave amplitude and timing is made by a second probe 54 placed on 
a fmger on the hand opposite that of the first probe. The pulse delay occurring between the 
two meeisurements is made. Alternatively, as shown in Figure 10, probes 52 and 54 can be 
placed on opposite temples of the patient to measure pulse wave values and delay. The 
probes can also be placed on the patient's ears. 

From this information alone, an estimate of pulse wave velocity at the aortic root 
could be made, by utilizing a table of normal values for the distance of the central 
anatomical difference. 

If a measurement of blood pressure is then made, one can perform the following 
calculation: 

Eq. 2 p = c*u*p 

where: 

c = pulse wave velocity; 
u = flow wave velocity; and 

p = the density of the blood (approximately 1 .055 grams/cm"'). 
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According to the invention, p and c are measured, and p is known. This allows one 
to solve for u, which is the flow wave velocity at the aortic root. This by itself is a 
measure of cardiac output. If one makes an estimate of aortic root diameter, one can then 
compute cardiac stroke volume. 

According to the invention, various well-known, conventional reconstruction 
techniques can be employed to convert or "transform" peripheral blood pressures and 
waveforms to the corresponding pressure and waveform at the aortic root. Ideally, the 
blood pressure at the aortic root should be used as the pressure term in Fung's equation. 

One can improve on the above determination in several ways. The first way is by 
additionally measuring the peripheral pulse wave velocity. To do this, measurement of 
pulse wave amplitude and timing is made by a first probe such as that shown in Figure 5. 
The probe is at a set known position relative to the heart. Another, simultaneous 
measurement of pulse wave amplitude and timing is made by a second probe placed on a 
finger on the hemd opposite that of the first probe. The pulse delay occurring between the 
two measurements is made. The respective peripheral pulse wave velocities are also 
computed. If the peripheral pulse wave velocities are different, it can be assumed that this 
is because of the different central anatomies from which the respective pulses traveled. 
This information alone may be enough to compute central pulse wave velocity from a 
table of normals. However, when combined with the pulse wave delay information, this 
data enables one to construct a function of pulse wave speed from the periphery back to 
the aortic root, thus giving another measure of central pulse wave velocity. 

Another method of the invention is to vary the position of the probes relative to the 
heart. If the first probe is at heeirt level and the second probe is raised above (with respect 
to the earth) heart level, the hydrostatic pressure of the blood vessels within the second 
probe will be lower than those within the first probe. In turn, in accordance with Fung's 
equation stated above, this means that the pulse wave velocity of the arterial vessels within 
the second probe will be lower them that in the arterial vessels within the first probe. This 
will change both the measured pulse delay between the two probes, and the measured 
peripheral pulse wave velocities. This creates additional measurements by which to 
compute central pulse wave velocity. 

According to the invention, changes in hydrostatic pressure are controlled by the 
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following equation: 

Eq. 3 p = p*g*h 

where: 

p = blood density; 

5 g — gravitational acceleration (approximately 9.8 meters/sec ); and 

h = height above a reference point (with respect to the earth). 
The difference in hydrostatic pressure between the vessels in two probes is thus 
governed completely by their difference in heights relative to the heart (referenced to the 
surface of the earth). Therefore, a known change in position produces a known change in 
10 hydrostatic pressure. 

According to the invention, the above measurements can be employed to derive a 
number of physiological properties. Preferably, the probes of the invention are connected 
to a controller to aid the data collection and analysis used to make the desired 
j.^ determination. The controller includes a computing device or standard personal computer 

t5 1 5 (PC) with a monitor. Included within the controller are algorithms for the calculation of 

tQ 

variables not measured directly. . 

^ For example, referring to Figure 1 1 , there is shown a circuit schematic for a one or 

5 two wavelength photo-plethysmograph. Emitters 56 and 58 and detector 60 are positioned 

l,i adjacent the tissue being measured, such as a finger 61 . Emitters 56 and 58 are driven by 

''3 20 drive circuitry 62, which is in turn governed by control signal circuitry 64. Detector 60 is 

□ connected to amplifier 66. The signal from amplifier 66 is sent to demodulator 68, which 

La 

is also synched to control signal circuitry 62. The signal from the demodulator 68 is sent 
to analog-digital converter 70. The desired computations are performed on the output 
from the converter 70 by signal processor 72 and the results sent to display 74. Emitters 

25 56 and 58 operate specific wavelengths, such as 805 nm, and may comprise light emitting 
diodes (LEDs) or laser diodes. Detector 60 preferably comprises a silicon photodiode. 
Such emitter-detector pairs are shown in Figures 2 and 3. 

Referring now to Figure 12, there is shown a schematic of an alternate embodiment 
of suitable circuitry. As with Figure 1 0, emitters 76 and 78 are connected via LED drive 

30 circuitry 79 and control signal circuitry 80 to demodulator 82. Signal from detector 84 is 
amplified at circuit block 86 and sent to demodulator 82. Output from demodulator 82 is 
sent to A/D converter 88. In addition, ECG leads 90 are connected to differential 
amplifier 92 and the signal is sent to converter 88. Output from converter 88 is processed 
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at block 94 and the results sent to display 96. A probe such as those shown in Figures 5 
and 6 may be used with the circuitry. The ECG leads are preferably silver/silver chloride 
or stainless steel. 

Yet another embodiment of the invention is shown in Figure 13. Emitters 98 and 
100 are connected via LED drive circuitry 101 and control signal circuitry 102 to 
demodulator 104. Signal from detector 106 is amplified at circuit block 108 and sent to 
demodulator 104. Output from demodulator 104 is sent to A/D converter 109. ECG leads 
1 10 are connected to differential amplifier 1 12 and the signal is sent to converter 109. 
Digit level sensor 114 and heart level sensor 116 are connected to amplifier 1 18 and the 
signal is sent to converter 109. Output from converter 109 is processed at block 120 and 
the results sent to display 122. 

Referring to Figure 14, there is shown a circuit schematic suitable for use with a 
probe having two physiczdly independent channels, such as the one shown in Figure 4. A 
first emitter-detector pair comprising emitters 124 and 126 and detector 128 are positioned 
adjacent the tissue being measured, such as a finger. A second pair comprising emitters 
132 and 134 and detector 136 are positioned a selected distance from the first pair. 
Emitters 124, 126, 132 and 134 are driven by drive circuitry 138, which is in turn 
governed by control signal circuitry 140. Signal from detector 128 is amplified by block 
142 and sent to demodulator 144. Independently, signal from detector 136 is amplified 
and demodulated at blocks 146 and 148, respectively. Output from demodulators 144 and 
148 is sent to sinalog-digital converter 150. The desired computations are performed on 
the output from the converter 150 by signal processor 152 and the results sent to display 



An alternative embodiment configured for use with a probe having two physically 
independent channels and an ECG lead, such as the one shown in Figure 7, is 
schematically illustrated in Figure 15. A^first emitter-detector pair comprising emitters 
156 and 158 and detector 160 are positioned adjacent the tissue being measured, such as a 
finger. A second pair comprising emitters 164 and 166 2ind detector 168 Eire positioned a 
selected distance from the first pair. Emitters 156, 158, 164 and 166 are driven by drive 
circuitry 1 70 which is in turn governed by control signal circuitry 1 72. Signal from 
detector 160 is amplified by block 174 and sent to demodulator 176. Independently, 
signal from detector 168 is amplified and demodulated at blocks 178 and 180, 
respectively. Output from demodulators 1 76 and 1 80 is sent to analog-digital converter 
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182. ECG leads 184 are connected to differential amplifier 186 and the signal is also sent 
to converter 1 82 The desired computations are performed on the output from the converter 
1 82 by signal processor 1 88 and the results sent to display 1 90. 

As one of ordinary skill in the art will appreciate, the placement of the various 
probes discussed above will affect the types of measurements that can be taken. As 
discussed above, Figures 9 and 10 show probes placed on opposite extremities to enable 
measurement of pulse wave delay. Figure 16 shows an embodiment of the invention with 
probe 52, such as in Figure 1 , placed on the digit, and a probe 54, such as in Figure 2, 
placed on the arm near the brachial artery. This could measure the pulse wave velocity in 
the arm (as well as pulse oximetry). A similar embodiment could measure pulse wave 
velocity in the leg. Figure 17 shows probes 52 and 54 placed on a finger and on a toe to 
measure the pulse wave delay. Figure 1 8 shows probes 52 and 54 placed on opposite 
digits and probe 55 placed on a toe. This allows measurement of the differential pulse 
wave delay between the fingers and toe, and allows calibration of the toe probe to be used 
in place of a finger probe (if only one finger probe could be used, such as in hand surgery). 
The use of appropriate probes also allows a diagnostic-quality ECG. Figures 19 and 20 
show probes 52 and 54 placed on opposite digits. One arm of the subject is placed at the 
level of the heart, while one arm is moved to different positions, both above and below the 
level of the heart. By generating different hydrostatic pressures in the vessels, the pulse 
velocity and hence pulse wave delay changes. In addition, the amplitude of the pulse 
wave, and amplitude of venous absorbance changes. This allows the additional 
computations of arterial blood pressure and venous pressure. Figure 21 shows probes 52 
and 54 placed on opposite digits and probe 55 placed on a toe. The differential hydrostatic 
pressures in the vessels allow measurements of pulse wave velocity and pulse wave delay, 
as well as arterial blood pressure and venous pressure. Use of probes with suitable ECG 
leads will also allow the invention to perform a diagnostic-quality ECG. In addition, heart 
rate and respiratory rate can be calculated, and cardiac output and several other 
cardiovascular characteristics computed. 

As discussed above, the controllers of the invention preferably output the results of 
the measurements and computations to a display. Referring now to Figure 22, there are 
shown two oscilloscope tracings fi"om pulse oximeter probes, such as those shown in 
Figure 1 , placed on the index fingers of both hands. The pulse wave delay is visable as the 
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slight phase difference between the two tracings. As the probes are at the same level, the 
pulse amplitudes are essentially identical. 

Referring to Figure 23, there are shown oscilloscope tracings for a probe on a hand 
placed at a level higher than the heart (Curve A) and a probe on a hand placed at a level 
lower than the heart (Curve B). The induction of a pressure differential between the two 
probes effects a change in the pulse delay. The change in pressure also correspondingly 
alters the pulse amplitudes. 

Referring now to Figure 24, there are shown oscilloscope tracings for a probe on a 
hand placed at a level lower than the heart (Curve A) and a probe on a hand placed at a 
level higher than the heart (Curve B). Here, the pulse delay has substantially reversed as 
have the pulse amplitudes. Figure 25 are oscilloscope tracings showing an 
electrocardiogram (Curve A) in conjunction with a pulse waveform (Curve B). 

The algorithms outlined below serve as examples, but modifications are possible to 
arrive at the indicated results, and are meant to be included within the spirit of this 
application. Various additional components of the device will be discussed in more detail 
below with reference to the following exemplary determinations. 

(D-1). Determination of Arterial Blood Pressure 

A probe, such as that shown in Figure 1 , is placed on an extremity that is moved in 
relation to the heart. As mentioned above, the hydrostatic pressure within the arteries and 
arterioles changes as a function of height with respect to the heart. Because of this, both 
the pulse wave velocity and pulse wave amplitude chemge as a function of probe height. 
These two parameters can be mapped against known distance above or below the heart. In 
this way, function curves of pressure vs. pulse wave amplitude and pressure vs. pulse 
wave velocity can be drawn. For example, a full excursion of the arm in a standing adult 
produces hydrostatic changes of greater than 50 cm of water in both directions. Using an 
arm and a leg, a gradient of well over 200 cm of water can be generated. This is a 
significant portion of the normal blood pressure range, and certainly enough to produce 
the function curves mentioned above. 

There is a huge amount of medical literature describing arterial behavior, so the 
curves can be extrapolated if necessary. These curves serve as calibration. 

It can thus be determined if "recalibration" is necessary ~ if either pulse amplitude 
or pulse wave velocity changes, and the other parameter does not change correspondingly. 
In other words, a shift on one curve should be matched by a corresponding shift on the 
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Other curve. If this shift does not occur as predicted, recahbration is required. Of course, 
the process of recalibration is the simple procedure outhned above. 

In a preferred embodiment, a first probe having a position sensor is placed level 
with the patient's heart. A second probe, such as one shown in Figure 8, having a position 
sensor and a pulse detector is placed on the patient's finger. The patient's arm is held out 
level with the heart so there is zero displacement between probes. Pulse amplitude is 
recorded from probe. The patient's arm is slowly raised, while pulse amplitude and 
relative displacement of probe are recorded. The hydrostatic pressure difference between 
probes is also computed. By comparing the recorded pulse £implitude to the hydrostatic 
pressure difference, a mathematiced function relating pressure to pulse amplitude can be 
derived. Preferably, circuitry similar to that shown in Figure 1 3 is used to aid the process. 
This process is repeated while lowering the arm back to heart level, then lowering the aim 
to below heart level and, finally, raising the arm back to heart level. Similar steps can be 
applied to measure pulse delay, pulse velocity and pulse contour. Also preferably, 
frequency or Fourier analysis may be used to make the requisite determinations. 

(D-2). Determination of Cardiac Outp ut 

Cardiac output can be determined by measuring delays in pulse arrival times in 
coupled organs or members on opposite sides of the body. In a preferred embodiment of 
the invention, probes such as those shown in Figure 1 , having sensors for detecting a 
patient's pulse are placed on opposite fingers of the patient. The patient positions both 
arms straight out from the side. The blood pressure of the patient can be determined either 
through conventional means or by the methods of the invention. The pulse delay between 
the two probes can be measured utilizing circuitry such as that shown in Figures 14 or 15, 
for example. The dicrotic notch of the pulse may be determined by standard methods, and 
used to calculate the ejection time based on the timing. The size of the aortic root can be 
estimated by standard means and the consequently the pulse distance differential at the 
aortic root. This allows the calculation of the pulse velocity c at the aortic route by the 
following equation: 

Eq. 4 c = (pulse distance)/(pulse delay) 

The value of c can then be used to determine the flow wave velocity based on the 
following equation: 

Eq. 5 p = c*u*p 
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where: 



c = pulse wave velocity; 
u = flow wave velocity; and 

p = density of the blood (approximately 1 .055 grams/cm^). 

According to the invention, cardiac stroke volume can be determined by 
multiplying the aortic root area by the flow wave velocity and by the cardiac ejection time. 
Cardiac minute output can be calculated by multiplying the cardiac stroke volume by the 
pulse rate. These steps can be augmented by raising and lowering the patient's arms with 
respect to each other to vary the pressure and the pulse wave velocity. 

Experiments by the inventor using a "single breath of oxygen" technique outlined 
below have shown that having the arms at different levels (and thus, different arterial 
blood pressures) results in different blood flows to the arms. The blood flow is thus 
obeying the classic fluidic version of Ohm's law: 



p = pressure 

f =flow 

r = resistance 

Because of this, varying the pressure by small amounts (large amounts would 
cause changes in autoregulation) can be used to calculate blood flow and thus cardiac 
output. Since the pulse arrives at the two hands at different times, if the "faster" arm is 
raised, the pulse wave velocity slows. The pressure difference necessary to make the 
pulse arrive in both hands simultaneously can be calculated. In view of the pressure-flow 
relationship above, the ratio of pressure should be the same as the ratio of flow. The time 
difference, or number of heartbeats, of the flow arrivals in the two hands when they are 
level with the heart can be computed. In turn, this gives the cardiac index, as the closer 
the flow arrival times are, the higher the cardiac output. Likewise, the greater the flow 
arrival time difference is, the lower the cardiac output. 

Alternatively, czu-diac output can be determined by placing probes such as those 
shown in Figure 5 on a patient's finger and toe. The probes measure oxygen saturation at 
each pulse. The oxygen saturation for each pulse at the first probe is compared to the 
oxygen saturation of that pulse and subsequent pulses at the second probe. With 
continuous monitoring, this allows the determination matching oxygen saturation, within 



Eq.6 



p = f^r 



where: 
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given tolerance limits, of the pulses from the probes. The patient's blood volume and the 
physical separation of the probes can be determined by standard methods. This allows the 
computation of cardiac stroke volume by dividing the blood volume displaced by the 
number of pulses. Then, the cardiac minute output can be calculated by multiplying the 
5 cardiac stroke volume by pulse rate. Circuitry such as that shown in Figures 11 or 12 is 
suitable for use with this embodiment. 

In yet another embodiment, the principles used above may be used to calculate 
cardiac output on the basis of flow wave propagation as opposed to pulse wave. As 
described above, the propagation of the flow wave (physical passage of blood through the 
10 arterial system) happens much more slowly than that of the pulse wave (celerity). Thus, 
blood that is ejected during a single cardiac cycle will arrive at the two hands at a slightly 
different time, while the time difference between the eurival at a hand versus a foot may be 
much longer, in the order of several seconds, 
p The difficulty lies in determining the blood that was ejected during the same 

P 15 cardiac cycle. Two parameters that can be used to determine this are the hemoglobin 
^f, oxygen saturation and the pH (acid-base balance). As will be discussed below in their 

respective sections, hemoglobin concentration and blood pH can be measured using the 

LH 

g probes and techniques of the invention. All blood that is ejected during a single cardiac 

cycle can be assumed to have the same saturation and pH, because blood is very well 

□ 20 mixed in the left atrium and left ventricle prior to ejection. Also, it can be assumed that the 

■1. J 

p saturation and pH change very little during the passage of the blood through the arterial 

M system, as the arteries behave like pipes, conveying the liquid but consuming or releasing 

minimal metabolites. 

Both the saturation and pH change slightly in a continuous fashion due to the 
25 respiratory cycle and other factors. In addition, large changes can be caused by maneuvers 
such as voluntary hyperventilation or breathholding. Because of this, a pattern in either the 
saturation or the pH can be recognized, allowing one to "match" the pulses which reflect 
the blood ejected during the same heartbeat. The difference in arrival times can then be 
used to compute the cardiac output. By way of illustration, referring to Figure 26, the first 
30 pulse train (Curve A) corresponds to the hand. The second pulse train (Curve B), which 
arrives later, corresponds to the foot. 

Figure 27 similarly illustrates this phenomenon for pulses of differing pH. The first 
pulse train (Curve A), which arrives earlier, corresponds to the hand. The second pulse 
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train (Curve B) corresponds to the foot. According to the invention, additional, known 
parameters in addition to oxygen saturation and pH can also be employed. 

As is well known in the art, oxygen saturation varies a small amount normally, even 
in the "steady state". The idea is that this variation can be used to "match" pulses that are 
of the same saturation, but are "out of phase." That is, they arrive at different parts of the 
body (opposite hands or hand and foot, for example) at different times. Changes in pH 
follow similar patterns. Referring to figures 28-30, there are shown induced changes in 
oxygen saturation caused by a single breath of oxygen, as measured by probes on the toe 
(Figure 28), left hand (Figure 29) and right hand (Figure 30), respectively. 

The concept of induced change in either oxygen saturation of pH is very important. 
Taking one of more breaths of oxygen will raise one's oxygen saturation. If one is already 
receiving oxygen therapy, then taking one of more breaths of a higher inspired 
concentration of oxygen ("richer mixture") will do. Conversely, breath holding, or 
breathing a gas other that oxygen, such as nitrogen, will decrease your oxygen saturation. 
For pH, breath holding or breathing a gas mixture containing carbon dioxide will lower 
your blood pH, whereas rapid breathing will raise your pH. As shown in Figures 28-30, 
the rise starts a few heartbeats earlier on the right hand compared to the left, and the foot is 
many beats later. 

For computing a value for cardiac output, the method is that an index can be 
generated based on patient's height or body surface area (BSA). The delay in the arrival 
of a given pulse (identified by saturation or pH) between different parts of the body 
generates an index for each pair of points. For exeimple, the longer the delay between 
hand and foot, the lower the cardiac output. The shorter the delay, the higher the cardiac 
output. 

With the addition of the induced saturation change (or pH), one can calculate all 
these indices more easily (since the change is more noticeable and easier to measure). 
Importzmtly, however, there is another index - the time between the breath inspiring the 
increased oxygen and the first heart beat in which the change is noted. This is the same 
principle - the shorter this time period, the greater the cardiac output, and the longer the 
time, the lower the flow. 

In summary, giving a subject a single breath of oxygen and continuously 
measuring the oxygen saturation in an extremity can be used for 3 measurements. First, 
the time or number of heartbeats from the inspiration to the first change of saturation in 
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the extremity (first pulse that has different saturation from baseline) can be used to 
compute the cardiac index (and thus cardiac output and stroke volume). Second, the time 
or number of heartbeats from the beginning of the change in saturation to the maximum 
change can be used to compute the cardiac ejection fraction. Finally, the time intervals of 
total time from inspiration until the saturation has returned to baseline, total time from 
beginning of saturation change to return to baseline and "rebound" in saturation after 
return to baseline can be used to compute blood volume. 

(D-3). Determination of Venous Saturation and Pressure 

Determination of arterial oxygen saturation can be determined by pulse oximetry 
and techniques well delineated in both the patent and medical literature. Hydrostatic 
changes as described in this application allow the determination of venous saturation and 
pressure as well. 

Place a probe, such as that shown in Figure 1 , on a fmger. Make measurements of 
both total absorbance and pulsatile absorbance. Raise the probe a known distance. Again 
measure both total absorbance and pulsatile absorbance. Both will be decreased. This is 
because the pulse amplitude is less because the arterial blood pressure within the probe is 
less (due to decrease in hydrostatic pressure). However, the total absorbance will also 
decrease, as the distending pressure in the venous system is less, and hence the veins and 
venules are smaller. All changes in absorbance can be assumed to be due to changes in 
blood volume. Saturation is calculated using the ratios of absorbance of distinct 
wavelengths. 

In one embodiment, the central venous pressure (CVP) can be estimated. A probe 
containing a position sensor is place level with a patient's heart. A second probe, such as 
the one shown in Figure 8, also comprising a position sensor is placed on the patient's 
fmger. The patient positions the arm so that the second probe is initially lower than the 
first probe. The total absorbance measured at the second probe is continuously monitored. 
The patient's arm is slowly raised, and the rate of change of absorbance of the second 
probe is computed with respect to the relative displacement to the first probe. When the 
rate of change changes by a predetermined amount representing an abrupt decrease, the 
arm position corresponding to the point of central venous drainage has been reached. The 
CVP can then be calculated by computing the hydrostatic pressure difference between the 
first probe and the second probe at that arm position. The circuitry shown in Figure 1 3 is 
suitable for use with this embodiment. 



29 



30 




ATTY DOCKET NO: \^^S-004C 

(D-4). Determination of Hemoglobin Concentration 
Pulse oximeters in clinical use measure hemoglobin absorbance at 2 different 
wavelengths, commonly 660 and 940 nanometers (nm). Ratiometric analysis of the 2 
absorbances at the 2 wavelengths allow determination of the relative amounts of the 2 
5 species of hemoglobin, oxyhemoglobin and deoxyhemoglobin, which normally constitute 
blood. Standard 2 wavelength oximeters assume that there only these 2 species. Though 
others may be present, they are usually so in very low concentrations. 

As discussed above, the basic equation governing the absorbance of light is Beer's 
Law: 

10 Eq. 7 A = E*C*L 

where: 

A = absorbance; 
E = the extinction coefficient; 
hi L = the path length of light travel; and 

□ 15 C = concentration of hemoglobin. 

to 

py The extinction coefficient of the medium through which the light is being 

transmitted can be thought of as transparency or opacity, 
g Using this nomenclature, the absorbance of the 660 nm wavelength is A^eo, and the 

absorbance of the 940 nm wavelength is A940. Also, oxy will refer to oxyhemoglobin and 

□ 20 deox to deoxyhemoglobin, 

p Writing out the absorbances in terms of Beer's Law provides: 

^'^ Eq. 8 A66O = E660oxy*Coxy*L+E660deox*Cdcox*L 

Eq. 9 A940 = E940oxy*Coxy*L+E940deox*Cdeox*L 

As sampling can be performed very quickly (optically and electronically), it can be 
25 assumed that path length L does not change between 2 successive measurements. Thus, L 
will cancel out when a division is performed. 

A66O E660oxy*Coxy*L+E660dcox*Cdeox*L 
A940 E940oxy * Coxy * L+E940deox*Cdeox* L 



E660oxy*Cpxy+E660deox*Cdeox 
E940oxy * Coxy'^E940deox * Cdeox 
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This can then be combined: 

E660oxy*Coxy*A940+E660deox*Cdeox*A940 = E940oxy*Coxy*A660 + E940dcox*Cdcox* A660 

With rearrangement: 

E660oxy*Coxy*A940 - E940oxy*Coxy*A660 — E940deox*Cdeox*A660 " E660deox*Cdeox* A940 
5 Coxy*(E660oxy*A940 - E940oxy*A66o) — Cdeox*(E940deox*A660 " E660dcox*A94o) 

Coxy/Cdeox = (E940deox* A660 - EeeOdeox* A940) / (E660oxy* A940 " E940oxy* A660) 

As the absorbances are measured and the extinction coefficients are known 
experimentally, this allows computation of Coxy /Cdeox- Designating this ratio Rs, the 
oxygen saturation of the blood (S) can be computed and the percentage of the blood that is 
10 oxyhemoglobin from the following: 
Eq. 10 S = Rs/(l+Rs) 

Unfortunately, Eq. 1 0 only allows determination of the ratio of the 2 species of 
1=.^ hemoglobin, not the total concentration of hemoglobin, more commonly termed the "blood 

^2 count". 

15 The inability of previous attempts to determine hemoglobin concentration from 

m 

absorbance measurements is because Beer's Law has 2 unknowns, C and L, with only a 
single measurement, A. As C and L are products, it is impossible to distinguish between 
3 them. Forming another equation by making another measurement of absorbance at, for 

example, a different path length, the new absorbance is simply a multiple of the first, and 

P 20 no new information is gained. 

H 

□ It will thus be appreciated that what is needed is the substitution for either C or L 

in a way that they are no longer linearly coupled. Since a non-linear result is desirable, 
multiplication or logarithms would be good choices to memipulate the equation. 
An absorbance measurement for hemoglobin can be expressed as: 
25 Eq. 1 1 A = C*L*k 

where: 

k = a value based on the extinction coefficients for the wavelength and the oxygen 
saturation of the blood. 

If there is only one species of hemoglobin present, k is the extinction coefficient 
30 for that wavelength. Dividing out this known value to provides: 
Eq. 12 A, =C*L, 



31 



ATTY DOCKET ]^™T-02-004C 



The absorbance measurements can then be taken at a different path length, 
resuhing in different values for A i as well as Li. However, it is not necessary to repeat 
measurements, as once absorbance values for a given path length are known, so are the 
values for all possible path lengths. This also holds for saturation. 

Eq. 13 A2 = C*L2 

Multiplying these Equations 1 1 and 12 provides: 

Ai*A2 = C*Li*C*L2 = C^*L,*L2 
Further, taking the logarithm of the new equation provides: 

Eq. 14 . log(A,*A2) = 2 log(C)+log(Li)+log(L2) 
Since the arithmetic operations may be performed in reverse order, this can be rewritten 
as: 

Eq. 1 5 log(L,)H-log(L2) = log(Ai * A2)-2 log(C) 

Taking logarithms first: 

Eq. 16 log(AO - log(C)+log(L,) 

Eq. 1 7 log(A2) = log(C)+log(L2) 

Next, multiplying these terms provides: 

Eq. 18 log(A,)*log(A2) = [log(C)+log(L,)] * [log(C)+log(L2)] 

Eq. 19 =log'(C)*log(C)*[log(L,)+log(L2)]+log(L,)*log(L2) 
Substituting the term [log(Li)+Iog(L2)] in Equation H2 provides: 

Eq. 20 log(A,)*log(A2) = log^(C)+log(C)*[log(A,*A2)-2 log(C)] 

+ log(L,)*log(L2) 

Rearrangement provides: 

Eq. 21 log(A,)*log(A2) = log(L,)*log(L2)+log(C)*log(A,*A2) 

- Iog(C)*2 log(C)+log'(L) 
= log(L,)*log(L2)+log(C)*log(A,*A2)-2 log^(C)+log^(C) 
= log(L,)*log(L2)+ log(C)*log(A|*A2)-log'(C) 
As above, taking the absorbance measurements at different path lengths, will give 
different values for Ai and A2, as well as Lj and L2. These new values are all multiples of 
our original A| and L|. Subscripts 3 and 4 are used to denote these new equations 
(measurements). The equations thus become: 

Eq. 22 log(A,)*log(A2) = log(L,)*log(L2)+log(C)*log(A,*A2)-log^(C) 

Eq. 23 log(A3)*log(A4) = log(L3)*log(L4)+log(C)*log(A3*A4)-log^(C) 
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Since the term log^(C) occurs in both equations, the term will be eliminated if the 
equations are subtracted from each other. Thus subtracting Equation 23 from Equation 22 
provides: 

Eq. 24 log(Ai)*log(A2)-log(A3)*log(A4) = log(L,)*log(L2) 

+ log(C)*log(A, * A2)-log\C)-[log(L3)*log(L4) 

+ log(C)*log(A3*A4)-log'(C)] 

= log(L,)*log(L2)-log(L3)*log(L4)+log(C)*log(Ai*A2) 

- log(C)*log(A3*A4) 

= log(L,)*log(L2)-log(L3)*log(L4)+log(C)*[log(A,*A2) 

- log(A3*A4)] 

These manipulations have broken the linear coupling of C and L, as L occurs as a 
square in this expression, while C has an exponent of 1 . Further, the A terms are in 
different combinations on opposite sides of the equation. In addition, multiplying either C 
or L by a given factor results in a different equation, not just a multiple of the same 
equation. New equations can thus be formed by taking multiples of either L or C. 
This equation also answers the questions: How can you tell C and L apart? How do you 
know that C is not half as large and L twice as large? The absorbance measurement is still 
the same. 

If LI is doubled and C halved, the A values should not change. Thus the left hand 
side of Equation 24 will not change. For the right hand side, halving C causes the value of 
the C term to decreeise by: 

Eq. 25 log 2* [log{A, * A2)-log(A3* A4)] 

Since the value of all the L will increase by log 2, the quantity of the L term will be: 

Eq. 26 log(L,*2)*log(L2*2)-log(L3*2)*log(L4*2) = [log(Li)+log 2] 

* [log(L2)+log2]-[log(L3)+log2]*[log(L4)+log2] 

= log(Li)*log(L2)+log2*(Iog L,+log L2)+log'2-[log(L3) 

* log(L4)+log2*(log L3+log L4)+log'2] 

= log(Li)*log(L2)-log(L3)*log(L4)+log2*(log L+log Lj-log La-log L4) 
Subtracting the original L term causes a change of: 

Eq. 27 log 2*(Iog L,+log Lz-log La-log L4) 

This is not equal to the change in the C term (if it were, L could be solved in terms 
of A), so C and L are not interchangeable. 
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Referring to Equation 24, it is noted that the log(C) term could be eUminated if 
new equations with 4 different path lengths were formed and the process repeated. There 
is nothing preventing this, as any number of equations at different path lengths can be 
formed, and all will be different. 
5 Eq. 28 log(A5)*log(A6)-log(A7)*log(A8) = log(L5)*log(L6)-log(L7) 

*log(L8)+log(C)*[log(A5*A6)-log(A7*A8)] 
To subtract Equation 28 from Equation 24, the coefficients of the log(C) terms 
must be equal, so Equation 28 is multiplied by the following term: 

[log(Ai*A2)-log(A3*A4)]/[log(A5*A6)-log(A7*A8)] 
10 To simplify, constants values that are known or already measured or determined 

are represented. 

Eq. 29 k,234 = log(Ai)*log(A2)-log(A3)*log(A4) 

Eq. 30 k5678 = log(A5)*log(A6)-log(A7)*log(A8) 

g Eq. 31 kd = [Iog(A,*A2)-log(A3* A4)]/[log(A5*A6)-log(A7* As)] 

□ 15 A new equation is formed by multiplying Equation 28 by kd. This ensures that the 

Co 

|»y log(C) terms are equal. Subtracting this new equation from Equation 24 eliminates the 

log(C) term: 

k,234-(kd*k5678) = log(Li)*log(L2)-log(L3)*log(L4) 

[: - kd*[log(L5)*log(L6)-log(L7)*log(L8)] 

□ 20 At this point, an equation with a constant (known number value) on one side, and 4 
p terms in log(L) on the other is obtained. Equations at 8 path lengths have been used, but 

these path lengths are all known multiples of each other. Thus, the equations can be 
expressed as: 

Eq. 32 L2 = Li*A2/A, 

25 Eq. 33 L3-Li*A3/Ai 

Eq. 34 L4 = L,*A4/Ai 

Eq. 35 L5 = L,*A5/Ai 

Eq. 36 L6 = Li*A6/Ai 

Eq. 37 L7 = L,*A7/A, 

30 Eq. 38 L8 = L|*A8/A| 
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The equation now becomes: 

Eq. 39 k,234-(kd*k5678) = log(L0*log(L,*A2/A0-log(L,*A3/A,) 

* log (L,*A4/Ai>kd*log(L,*A5/Ai)*log(L,*A6/Ai)+kd 

* log(L,*A7/Ai)*log(L,*A8/Ai) 

5 =log(L,)*[log(L,)+log(A2/A0]-[log(L,)+log(A3/A,)]*[log(L,) 

+ log(A4/A,)]-kd*[log(L0+log(A5/A,)]*[log(L,)+log(A6/A0] 
+ kd*[log(L,)+log(A7/AO]*[log(L,)+log(A8/A,)] 
= log'(L,)+log(L,)*log(A2/A,)]-[log' (L,)+[log(L,)*(log(A3/Ai) 
+ log(A4/A,)]+log(A3/A,)*log(A4/Ai)]-kd*[log' (Li)+log(Li) 
10 * (log(A5/Ai)+log(A6/Ai)]+log(A5/A,)*log(A6/Ai)]+kd*[log' (Li) 

+ [log(Li)*(log(A7/Ai)+log(A8/Ai)]+log(A7/A,)*log(A8/A,)] 
As noted, the log^ terms have coefficient 1 . Further, since 2 terms are added while 
2 are subtracted they will be eliminated. 

Combining all the A terms (which are known) that occur as coefficients, creates 
C3 15 kAc- Then, combining all the A terms that occur alone (kAa)- the equation becomes: 

fU Eq. 40 ki234-(kd*k5678) = log(Li)*kAc+kAa 

1X3 

where: 

5 log(Li) = [ki234-(kd*k5678)-kAa]/kAc 

^s. Li - antilog{[ki234-(kd*k5678)-kAa]/kAc} 

P 20 This leaves a single unknown, Li. Note that all A values are different, so neither side 

p of the equation should be zero. If this does occur, choice of one or more different A 

values will guarantee a non-zero result. 

Since the equation has been solved for Li, returning to the original equation provides: 
Eq.41 Ai = C*Li 

25 Eq. 42 C = A,/Li 

Now, the concentration of hemoglobin, C, can be solved. An example of a suitable 
computational algorithm is provided in Figure 33. 

A similar analysis can be performed, varying C instead of L. Maintaining the 
assumption of only 2 species of hemoglobin, then the wavelength absorbances can be 
30 expressed as: 

Eq. 43 A660 = A660oxy+A660deox 

Eq. 44 A940 = A940oxy+A940deox 
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Furthermore, 

Eq. 45 Aoxy = A660oxy+A940oxy 

Eq. 46 Adcox = A660deox+A940deox 



and 



Eq. 47 Aoxy = E660oxy*Coxy*L+E940oxy*Coxy*L 

Eq. 48 Adcox = E660deox*Cdeox*L+E940deox*Cdeox*L 



Thus: 



Eq. 49 A660oxy E660oxy*Coxy*L 

10 A660dcox E660deox*Cdeox*L 

Since the ratios of Coxy-Cdeox are known, the equations can be written (cancehng L 
in the process): 

Eq. 50 A660oxy E660oxy*Ctot*S E660oxy*S 

p A660deox E660deox*Ctot*(l-S) E660deox*(l-S) 

ci 

Knowing saturation S and the extinction coefficients, permits the ratio, defined as 
rU R660, to be solved. Returning to the absorbance equations Aeeooxy and A^eodcox^ 

En 20 A660oxy = A660*[R66o/(l+R66o)] 

A660deox = A660-A660oxy 

M The process can be repeated for the 940 nm absorbance, A940. In summary, each of 

^.^ the individual absorbEmces, i.e., Aeeooxy, Aeeodeox, A940oxy, A940deox has been determined. 

^3 Combining them as above provides: Aoxy, Adeox, which can be rewritten: 

25 Eq. 51 Aoxy = Coxy* L*(E660oxy"^E940oxy) 

Eq. 52 Adeox ~ Cdeox*L*(E660deox+E940deox) 

At this point it is convenient to divide out the known extinction coefficient values 
to get A'oxy and A'deox, which remain known values: 

Eq. 53 A oxy ~ Coxy*L 

30 Eq. 54 A'dcox = Cdeox*L 

As above, these equations can be manipulated with multiplication and logarithms. 
Multiplying first: 

Eq. 55 A oxy*A deox ~ Coxy *L*Cdeox*E — Coxy*Cdeox*L' 

Then taking the logarithm: 

35 Eq. 56 log(A'oxy*A'deox) = log(Coxy)+log(Cdeox)+2 log(L) 
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This can be rewritten as: 

Eq. 57 log(Coxy)+log(Cdcox) = log(A'oxy* A'deox)-2 log(L) 

The operations can now be performed in reverse order. 
Taking the logarithms first: 

Eq. 58 log(A'oxy) = log(Coxy)+log(L) 

Eq. 59 log(A'deox) = log(Cdeox)+log(L) 

Now multiplying these terms: 

Eq. 60 log(A'oxy)*log(A'deox) = [log(Coxy)+log(L)]*[iog(Cdeox)+log(L)] 

= log(Coxy)*log(Cdeox)+log(L)* [log(Coxy)+log(Cdeox)]+log'(L) 

Now substituting for log(Coxy) + log(Cdeox) from Equation 57 into Equation 60: 

Eq. 61 log(A'oxy)*log(A'deox) = log(Coxy)*log(Cdeox)+log(L) 

*[log(A'oxy*A'deox)-2 10g(L)]-f log2(L) 

Rearranging this equation provides: 

Eq. 62 log(A'oxy)*log(A'dcox) = log(Coxy)*log(Cdcox)+log(L)*log(A'oxy 

* A'deox)-log(L)*2 log(L)+log2(L) 

= log(Coxy)*log(Cdeox)+log(L)*log(A'oxy* A'deox )-2 log\L)+log\L) 
= log(Coxy)*log(Cdeox)+log(L)*log(A'oxy*A'deox)-log^ (L) 

Taking absorbance measurements at a different saturation gives different values for 
A'oxy and A'deox, as well as Coxy and Cdeox* As above, once absorbance values for a given 
saturation are known, they are known for all possible saturations. 

Therefore, A'oxy and A'deox can be calculated for a different saturation, and 
performing the multiplication and logarithm manipulations will generate a new equation. 
The original values of A'oxy and A'deox are specified as Aoxyi and Adeoxi, and the new 
values, Aoxy2 and Adeox2- 

The equations thus become: 

Eq. 63 log(Aoxyl)*log(Adeoxl) log(Coxyl)*log(Cdeoxl)+log(L) 

* l0g( Aoxy 1 * Adeox I log^(L) 

Eq. 64 log(Aoxy2)*log(Adeox2) = log(Coxy2)*log(Cdeox2)+log(L) 

* log(Aoxy2*Adeox2)-log^(L) 

Since the term log^(L) occurs in both equations, it will be eliminated by subtracting 
the equations from each other. Thus, subtracting Equation 64 from Equation 63 provides: 
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Eq. 65 l0g(Aoxyl)*l0g(Adeoxl)-l0g(Aoxy2)*l0g(Adeox2) = log(Coxyl) 

* log(Cdeoxl )+log(L)* 10g(Aoxy 1 * Adeox 1 )-log^(L)-[log(Coxy2) 
*log(Cdeox2)+log(L)*log(Aoxy2*Adcox2)-log^(L) 

= log(Coxyl)*10g(Cdcoxl)-log(Coxy2)*log(Cdcox2)+log(L) 

* log(Aoxyl*Adeoxl)-log(L)*log(Aoxy2*Adcox2) 

Eq. 66 = 10g(Coxyl)*log(Cdeoxl)-log(Coxy2)*10g(Cdeox2)+10g(L) 

* [log(Aoxyl*Adeoxl)-log(Aoxy2* Adeox2)] 

Referring to Eq. 66, it can be seen that the log(L) term could be eHminated 
according to the invention, by forming new equations with 2 different saturations and 
repeating the process. Again, any number of equations at different saturations can be 
formed, and all will be different. 

Using reference numerals 3 and 4 to denote these new relationships, provides: 

Eq. 67 log(Aoxy3)*log(Adeox3)-log(Aoxy4)*log(Adeox4) = log(Coxy3) 

* log(Cdeox3)-log(Coxy4)*log(Cdcox4)+log(L)*[log(Aoxy3*Adeox3 ) 

- log(Aoxy4*Adeox4)] 

To subtract Equation 67 from Equation 66, the coefficients of the log(L) terms 
must be equal, so Equation 67 can be multiplied by the following: 

[log(Aoxyl*Adcoxl)-log(Aoxy2*Adeox2)]/[log(Aoxy3*Ade(>x3)-log(Aoxy4*Adeox4)] 

To simplify, representing by values of constants that are known or already 
measured or determined: 

Eq. 68 ki2 = log(Aoxyl)*log(Adeoxl)-log(Aoxy2)*log(Adeox2) 

Eq. 69 k34 = log(Aoxy3)*log(Adeox3)-log(Aoxy4)*log(Adeox4) 

Eq. 70 = [log(AoxyI *Adeoxl)-log(Aoxy2* Adeox2)]/[log(Aoxy3* Adeox3) 

- log(Aoxy4*Adeox4)] 

A further equation is formed by multiplying Equation 67 by km. This ensures that 
the log(L) terms are equal. Subtracting this new equation from Equation 66 eliminates the 
log(L) term: 

Eq. 71 k,2-(k„,*k34) = log(Coxyl)*log(Cdeoxl)-log(Coxy2)*log(Cdeox2) 

-km*[log(Coxy3)*log(Cdeox3)-log(Coxy4)*10g{Cdeox4)] 

At this point, the equation has a constzint on one side, and 4 terms in log^C on the 
other. Equations at 4 values of saturation have been used, but these saturation values are 
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all known. Thus: 




Eq. 72 


^oxy 1 


^^lOl 1 


Eq. 73 






Eq. 74 






riq. I J 


<-deox2 - 


- Mot (l-2>2) 


Eq. 76 


Coxy3 ~ 


Ctot*S3 


Eq. 77 


Cdeox3 ~ 


= Ctot*(l-S3) 


Eq. 78 


Coxy4 ~ 


Ctot*S4 


Eq. 79 


Cdeox4 ~ 


= Ctot*(l-S4) 



10 The equation now becomes: 

Eq. 80 k,2-(k,„*k34) = log(C,o.*S,)*log(Cu„*(l-S,))-log(C,o,*S2) 

* log(Cu«*(l-S2)>- k„,*log(C,o.*S3)*log(CK,.*(l-S3))+k„ 

i.4 * l0g(Cu«*S4)*l0g(Cu„*(l-S4)) 

Inspecting the 4 terms on the right hand side of this equation individually: 

h 15 Eq. 81 log(C,„,*S,)*log(C«„*(l-Si)) = [log(C,„,)+log(S,)]*[log(Qo.) 

fy +iog(i-s,)] 

=log\C,o,)+iog(C,o,)*log(S,)+log(Cto,)*Iog(l-S,)]+log(S,)*log(l-S,) 
= log'(C,o,)+log(C,ot)*[log(Si)+log(l-Si)]+log(Si)*log(l-S,) 
[7 Note that the log^ term has 1 as coefficient. Since there are 4 symmetrical terms, 2 

C3 20 adding and 2 subtracting, the log^ terms will cancel out i.e., 
£3 Eq.82 k,Hk,n*k34) = log(Ctot)*[log(Si)+log(l-S,)]+log(S,) 

^' * lOg(l-Sl)-lOg(C,o0*[lOg(S2)+log(l-S2)]-log(S2) 

♦ log( 1 -S2)-k™* [log(C,o,)* [log(S3)+log( 1 -S3)]+log(S3) 

» l0g( 1 -S3)]+k„* [l0g(Cu,t)* [l0g(S4)+l0g( 1 -S4)]+l0g(S4) 
25 * log(l-S4)] 

All terms containing only log(S) or Iog(l - S) are known, so they can be collected 
and replaced with kj. Further, collecting all terms with C,ot provides: 

Eq. 83 k,2-(k„,*k34) = log(Cu,0*[[log(S,)+log(l-S,)]-[log(S2) 

+ log(l-S2)]-k„*[log(S3)+log(l-S3)]+k„*[log(S4) 

30 + log(l-S4)]]+ks 
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Again, all tenns with log(S) or log(l-S) are known, so they can be collected and 
replaced with kcs. Note that kcs cannot be zero, as the saturations are all different. The 
equation thus becomes: 

Eq. 84 ki2-(k„,*k34) = log(Ctot)*kcs+ks 

where: 

ki2 - (km*k34) - ks = log(Ciot)*kcs 

log(Ctot) = [k,2-(kn,*k34)-ks]/kcs 

Ctot = antilog{[ki2-(km*k34)-ks]/kcs} 

In summary, Ctot, the concentration of hemoglobin, has been determined based 
entirely on known or measured values. It will be appreciated that, while this analysis 
£issumes only 2 species of hemoglobin, the determinations can be made with additional 
species present by utilizing additional wavelengths of light. 

If the saturation is known, the determination can be made using a single absorbance 
measurement from a single wavelength of light. For illustration, absorbance employing a 
wavelength of 660 nm is detailed below. 

Eq. 85 A660 = A660oxy+A660deox 

^660 ~ E660oxy*Coxy*L+E660deox*Cdcox*L 

AggOoxy E660oxy * Coxy * L 

A660deox EeeOdeox* Cdcox* L 

Rewriting this equation to incorporate the known saturation S gives: 

Eq. 86 A660oxy E660oxy*Ctot*S E660oxy*S 

AeeOdeox E660deox*Ctot*(l-S) E660deox*( 1 "S) 

As above, this equation can be solved, providing the ratio Reeo- Returning and 
solving for the actual absorbances Aseooxy and Aeeodeox provides: 

Eq. 87 A660oxy = A660*[R660/(1+R660)] 

A660deox ~ Aeso-Aeeooxy 
Finally, dividing out the known extinction coefficient values provides A'660oxy and 

A'660deox* 

Eq. 88 A'660oxy = Coxy*L 

Eq.89 A'660deox = Cdcox*L 
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At this point, the same algorithm as above will work for the single wavelength, 
allowing the substitution of the 660 nm absorbances for the total absorbances. That is, 
substituting A'^gooxy for A' oxyj etc. It will also be appreciated that performing ansdysis by 
both methods (manipulating C or manipulating L) independently can be used to verify 
results obtained by either technique. 

An altemative measurement method utilizing venous absorbance follows. The 
technique is equally suitable for arterial measurements. Two probes, such as shown in 
Figure 8, are placed on two separate digits. These probes should be of fixed or at least 
known dimensions. The probes can be on one hand or both hands, as in Figure 9. Raising 
the hand well above the heart causes the venous absorbance signal to fall to near zero, as 
height of the probes will cause the pressure in the venous system within the probes to be 
below the CVP. Slowly lowering the digits will cause the veins to no longer be collapsed 
at some point. While this should happen at the same point for both probes, it is not 
essential to the technique. As the lowering of the digits continues, the volume and 
pressure will rise above zero, and will continue rising as the digits are lowered. The 
important concept is that the absorbances will not be the same in both probes. Two 
dissimilar fingers, such as a thumb and fifth finger can be chosen, but in any case 
absorbances will not be identical even between similar fingers. Circuitry such as that 
shown in Figure 13, may be used to facilitate the measurements. 

Since the dimensions of the probes are known, the concentration C of the blood 
within the probes can be imputed. Of course, this C will be much higher than the true 
concentration since the space within the probes is occupied by other tissue and only a 
small part is blood. If H represents the hydrostatic distance that the probes have moved 
since the beginning of venous absorbance (this is convenient to establish a zero point for a 
graph), then a graph of H vs. C will be different for the two probes, as the amount of 
absorbing blood within the two probes is different. Given that the true concentration of 
hemoglobin in the blood is the same throughout the vasculeir system, there is a fiirther 
distance D (a fraction of H) that a first probe with the lower concentration would have to 
be physically lowered in order to achieve the same concentration as the second probe. 
Preferably, the dimensions of the probes should be minimized to increase the accuracy of 
the concentration estimates, since the path length of the probes is much greater than the 
path length of blood. For example, halving the diameter of the probes doubles the 
measured concentration. 
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Since the true path length of non- scattering blood absorbance is very small, the 
point where H and (H+D) are in the same ratio as the imputed concentrations yields the 
true concentrations. That is, if the concentration ratios are 1 :2 (twice as much absorbance 
in one probe), the true concentrations would be where H and D are equal. The actual 
concentration of hemoglobin is that at point (H,C). The further length D represents the 
additional non-blood path length through the other probe. Note that this C occurs at only 
one point. Also note that H is not the true path length, but C is in fact the true 
concentration. 

Advantageously, this technique can be repeated many times for one arm 
movement, can be used for both arterial and venous measurements, and arterial and 
venous measurements can be done simultaneously, and checked against each other. 

In another embodiment of the invention, hemoglobin concentration can be derived 
by taking ratios of absorbemce at different pressures. Pressure is equal to dynes divided by 
the square of the path length, i.e., 

Eq. 90 P = DO? 

and, as indicated above, Beer's Law is: 

Eq. 91 A = E*C*L 

Multiplying P by A^ provides: 

Eq. 92 E^*C^*D 

Successive measurements give E^*C^*Di, E^*C^*D2, . . . E^*C^*Dx from which 
ratios can be formed. Next, the values P + A^ and P"' + A^ are known quantities, and 
substituting Beer's Law and the pressure to length relationship allows the latter to be 
written as: 

Eq. 93 lW+E^*C^*l2 
Dividing by L^ provides: 

Eq. 94 L^*(D-Ve^*C^) 

As before, successive measurements gives values Li^*(Dr^+E^*C^), L2^*(D2" 
^+E^*C^), . . . Lx^*(Dx'^+E^*C^). Since the path length ratios are known, the values can be 
divided by L and the successive values can be divided to give the ratio: 

Eq. 95 (Di-^+E^*CV (D2"^+E^*C^) - m 

Further, Di'^/Dj'^ can be expressed in terms of P2/D1, which, as discussed above, is 
a known ratio. As E^*C^ is a constant, equation 95 is taking a known ratio, adding a 
constant to the numerator and denominator, and arriving at another known ratio. Thus, the 
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value of the constant can be determined uniquely and precisely. As E is known, this 
allows determination of C, the concentration of hemoglobin. Once C is determined, 
A = ECL allows computation of L, the path length and determination of L allows the 
computation of D as validation. 

(D-5). Determination of pH and Acid-Base Balance 

This probes and methods of this invention allow the pH and acid-base balance of a 
patient's blood to be determined noninvasively, continuously and in real time. Generally, 
a probe such as those described above, containing one or more sources of emitted light is 
used, for example as shown in Figures 5 and 8. Preferably, the invention utilizes light in 
the visible and near-infrared spectra, but other wavelengths could be successfully used. 
Detection of transmitted light is made by a photodetector. Suitable circuitry is shown in 
Figures 1 1 and 14, for example. Determination of absorbance of the emitted light is done 
using Beer's law. The emission, detection, and determination of absorbance described are 
all similar to that used in pulse oximeters, common in modem hospital settings. Using 
either pulse oximetry (filtering out all non-pulsatile signal components) or conventional 
spectrophotometry, identification of methemoglobin, oxyhemoglobin, deoxyhemoglobin, 
and carboxyhemoglobin, species of human hemoglobin, is made by utilizing known 
absorbance spectrums, such as those shown in Figure 34. 

The absorbance spectra of oxyhemoglobin, deoxyhemoglobin, and 
carboxyhemoglobin, the more common hemoglobin species, are conventionally 
understood to be essentially unaffected by acid-base balance over the most commonly 
used spectrophotometric range of 650 - 1000 nm. However, all proteins and organic 
molecules are affected by pH to some extent. Figure 35 shows the absorption spectrum of 
hemoglobin A (deoxyhemoglobin) between 500 and 650 nm for 4 different pH values - pH 
6, 7, 8, and 9. It can be clearly seen that pH does have a significant effect on the 
absorption spectrum. Thus, the spectrum of hemoglobin species is in fact altered by pH in 
the clinical range. The changes are small, and have largely been ignored by the prior art 
because they are too small to be considered clinically important for oxygen saturation 
measurements. 

While many wavelengths may be affected, specific alterations have been found at 
wavelengths 535, 577, and 600 nm. For example, oxyhemoglobin shows absorption peaks 
at 535 and 577 nm that are greater at higher pH, while the trough near 600 nm is slightly 
deeper as illustrated in Figure 36, The upper panel is a difference spectrum magnified 
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10- fold to more clearly illustrate the absorption pattern. The total hemoglobin 
concentration in both curves was normalized arbitrarily at a wavelength of 506 nm. It is 
estimated that this produces a difference in calculated saturation of about 3% from pH 7.6 
to pH 6.8. Thus, saturation mesisurements accurate to 0.1% should define pH to within 
0-05 pH units. 

Ratiometric techniques similar to pulse oximetry can be used to exploit these 
spectral differences in order to determine pH of the blood. As pH has relatively little 
effect at the near infrared wavelengths, performing pulse oximetry using two wavelengths 
such as 660 and 940 can be used to determine the "true" oxygen saturation of the blood 
(for now, assuming only the two species oxyhemoglobin and deoxy hemoglobin). An 
additional simultaneous measurement of at least one other wavelength ~ one known to be 
affected by pH, for example 535 or 577 nm — is performed. By making pair- wise 
comparisons of the ratios at all the wavelengths, the pH can be measured from the 
saturation values. Preferably, a multi-wavelength spectrophotometer allows multiple 
comparisons. 

Similarly, there are small pH-dependent changes in the absorption spectra of 
carboxyhemoglobin and fetal hemoglobin. A pH effect has been determined on the 
carboxyhemoglobin (COHb) absorption spectrum. There is a bias of about 3% COHb per 
pH unit if pH is not taken into account. Figure 37 illustrates an absorption spectrum and 
shows that the absorption peaks at 535 and 577 nm are higher, and the trough near 600 nm 
is slightly lower, at the higher pH value. Thus, the use of ratiometric determinations of 
predicted versus actual fraction of carboxyhemoglobin can be used to calculate pH of 
blood noninvasively using oximetric techniques. Differences have been found in the 
absorption spectra of species of fetal hemoglobin due to pH. especially fetal 
carboxyhemoglobin. This is illustrated in Figures 31 and 32. Accordingly, noninvasive 
pH determinations can be made from its absorption ratios. 

As mentioned above, the absorbance spectra of oxyhemoglobin, deoxyhemoglobin, 
and carboxyhemoglobin, the more common hemoglobin species, are affected relatively 
little by acid-base balance. Another embodiment of the invention exploits the optical 
spectrum of methemoglobin, which is very sensitive to pH in the physiologic range. 
Specifically, at pH 7 and below, the spectrum of aquomethemoglobin with maximum at 
wavelength 630 nm is observed. At alkaline pH, the spectrum of hydroxymethemoglobin 
with disappearance of the maximum at 630 nm is seen. In terms of electron paramagnetic 
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resonance spectra reflecting the electronic spin state of the iron atom, in 
aquomethemoglobin the iron atom has a "high spin" state and is strongly paramagnetic. 
With increasing pH and the formation of hydroxymethemoglobin, the iron atom goes to a 
low spin state, accompanied by a decrease in atomic radius. Since methemoglobin exists 
as an equilibrium mixture of structures, the relative amount of each can be used to 
determine the pH of the medium, blood. 

Referring now to Table I, there is shown a table by which the extinction 
coefficients for methemoglobin for selected wavelengths of light can be determined, since 
the absorbance varies with the acid-base balance of the blood. 

TABLE I 



Species 


Wavelength 

540 560 570 576 640 


Oxyhemoglobin 


1.53 0.906 1.23 1.65 0.015 


deoxyhemoglobin 


1.08 1.34 1.16 1.01 0.115 


carboxyhemoglobin 


1.44 1.22 1.51 1.13 0.023 






pH 


Methemoglobin 


6.2 


0.609 0.373 0.355 0.358 0.412 


6.4 


0.610 0.374 0.356 0.360 0.410 


6.6 


0.613 0.379 0.365 0,370 0.409 


6.8 


0.623 0.389 0.380 0.385 0.406 


7.0 


0.638 0.405 0.400 0.406 0.401 


7.2 


0.658 0.424 0.425 0.433 0.393 


7.4 


0.682 0.450 0.455 0.465 0.380 


7.6 


0.710 0.481 0.492 0.507 0.360 


7.8 


0.746 0.520 0.537 0.558 0.324 


8.0 


0.788 0.567 0.590 0.620 0.306 


8.2 


0.835 0.616 0.647 0.682 0.280 


8.4 


0.882 0.665 0.703 0.745 0.254 


8.6 


0.926 0.710 0.755 0.800 0.230 


8.8 


0.964 0.749 0.801 0.850 0.208 



From the information in Table I, Tables IIA and TIB can be created, which gives 
the ratios of the extinction coefficients for methemoglobin at the various pH values. To 

determine the pH using Tables IIA and IIB, one first determines either the absorbance due 
to methemoglobin at each wavelength of interest, or the percentage of methemoglobin 
compared to other hemoglobin species. These determinations can be done using 
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techniques described elsewhere in this application, or in other works. Once this has been 
done, the ratio of absorbances at the different wavelengths can be used to compute the pH. 

For example, if the ratio of absorbance of methemoglobin for a 540 nm wavelength 
light source to the absorbance of methemoglobin for a 560 wavelength light source is 
1.633, the pH is 6.2. If the ratio is instead 1.516, the pH is approximately 7.4. Similar 
calculations can be made for other wavelengths, and the results compared to ensure 
accuracy. 



TABLE IIA 







Methemoglobin 




pH 


540/560 


540/570 


540/576 


540/640 


560/570 


6.2 


1.633 


1.715 


1.701 


1.478 


1.051 


6.4 


1.631 


1.713 


1.694 


1.488 


1.051 


6.6 


1.617 


1.679 


1.657 


1.499 


1.038 


6.8 


1.602 


1.639 


1.618 


1.534 


1.024 


7.0 


1.575 


1.595 


1.571 


1.591 


1.013, 


7.2 


1.552 


1.548 


1.520 


1.674 


0.998 


7.4 


1.516 


1.499 


1.467 


1.795 


0.989 


7.6 


1.476 


1.443 


1.400 


1.972 


0.978 


7.8 


1.435 


1.389 


1.337 


2.302 


0.968 


8.0 


1.390 


1.336 


1.271 


2.575 


0.961 


8.2 


1.356 


1.291 


1.224 


2.982 


0.952 


8.4 


1.326 


1.255 


1.184 


3.472 


0.946 


8.6 


1.304 


1.226 


1.158 


4.026 


0.940 


8.8 


1.287 


1.203 


1.134 


4.635 


0.935 
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TABLE IIB 







Methemoglobin 




pH 


560/576 


560/540 


570/576 


570/640 


576/640 


6.2 


1.042 


0.905 


0.992 


0.862 


0.869 


6.4 


1.039 


0.912 


0.989 


0.868 


0.878 


6.6 


1.024 


0.927 


0.986 


0.892 


0.905 


O.o 


1 ni n 

1 .Ul u 


0.958 


0.987 


0.936 




7 0 




1.010 


0.985 


0.998 


1 019 




0 Q7Q 
\j.y / y 


1.079 


0.982 


1.081 


1 1 09 


1 4 




1.184 


0.978 


1.197 


1 994 


7.6 


0.949 


1.336 


0.970 


1.367 


1.408 


7.8 


0.932 


1.605 


0.962 


1.657 


1.722 


8.0 


0.915 


1.853 


0.952 


1.928 


2.026 


8.2 


0.903 


2.200 


0.949 


2.311 


2.436 


8.4 


0.893 


2.618 


0.944 


2.768 


2.933 


8.6 


0.888 


3.087 


0.944 


3.283 


3.478 


8.8 


0.881 


3.601 


0.942 


3.851 


4.087 



In these embodiments of the invention, it is not necessary to compute the actual 
concentration of any hemoglobin species. Indeed, merely the ratio or relative amount 
compared to other species yields sufficient information to determine blood pH and acid- 
base balance. In general, the separate identification of n species will require at least n 
different wavelengths of light to be used. If desired, more tlian the minimum number of 
wavelengths may be used. 

These embodiments of the invention can be supplemented by considering the effect of 
temperature oil the pH of blood. Temperature not only affects the hemoglobin oxygen 
dissociation curve, but has an independent effect on pH. Figure 38 is a graph of the normal 
hemoglobin-oxygen dissociation curve (HODC). Figure 39 shows the effect of temperature 
and pH on the HODC of Figure 38. Thus, measurements of spectra can be taken at more that 
one temperature (and thus more than one pH), improving the accuracy of determinations. The 
Henderson-Hasselbach equation, describes the dissolution of an acid in terms of pH, pK 
(dissolution or dissociation constant), and the concentrations of the acid and its salt or base. 
The solubility, X, of carbon dioxide (CO2) is temperature-dependent and the pK for CO2 also 
depends on temperature. Thus, for CO2 the Henderson-Hasselbach equation becomes pH = 
pK + Iog([HC03"]/ XPCO2), or an alternate form can be used. The temperature dependence 
of the pK for carbonic acid is shown in Tables IIIA - HIE below. 
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# 



TABLE IIIA 


Temperature 


PH 

7.6 7.59 7.58 7.57 7.56 7.55 7.54 7.53 


40 


6.07 6.0705 6.071 6.0715 6.072 6.0725 6.073 6.0735 


39 


6.07 6.0705 6.071 6.0715 6.072 6.0725 6.073 6.0735 


38 


6.08 6.0805 6.081 6.0815 6.082 6.0825 6.083 6.0835 


37 


6.08 6.0805 6.081 6.0815 6.082 6.0825 6.083 6.0835 


36 


6.09 6.0905 6.091 6.0915 6.092 6.0925 6.093 6.0935 


35 


6.1 6.1005 6.101 6.1015 6.102 6.1025 6.103 6.1035 


33 


6.1 6.1005 6.101 6.1015 6.102 6.1025 6.103 6.1035 


30 


6.12 6.1205 6.121 6.1215 6.122 6.1225 6.123 6.1235 



TABLE IIIB 



Temperature 


PH 

7.52 7.51 7.5 7.49 7.48 7.47 7.46 7.45 7.44 


40 


6.074 6.0745 6.075 6.0755 6.076 6.0765 6.077 6.0775 6.078 


39 


6.074 6.0745 6.075 6.0755 6.076 6.0765 6.077 6.0775 6.078 


38 


6.084 6.0845 6.085 6.0855 6.086 6.0865 6.087 6.0875 6.088 


37 


6.084 6.0845 6.085 6.0855 6.086 6.0865 6.087 6.0875 6.088 


36 


6.094 6.0945 6.095 6.0955 6.096 6.0965 6.097 6.0975 6.098 


35 


6.104 6.1045 6.105 6.1055 6.106 6.1065 6.107 6.1075 6.108 


33 


6.104 6.1045 6.105 6.1055 6.106 6.1065 6.107 6.1075 6.108 


30 


6.124 6.1245 6.125 6.1255 6.126 6.1265 6.127 6,1275 6.128 


TABLE IIIC 


Temperature 


PH 

7.43 7.42 7.41 7.4 7.39 7.38 7.37 7.36 7.35 


40 


6.0785 6.079 6.0795 6.08 6.0805 6.081 6.0815 6.082 6.0825 


39 


6.0785 6.079 6.0795 6.08 6.0805 6.081 6.0815 6.082 6.0825 


38 


6.0885 6.089 6.0895 6.09 6.0905 6.091 6.0915 6.092 6.0925 


37 


6.0885 6.089 6.0895 6.09 6.0905 6.091 6.0915 6.092 6.0925 


36 


6.0985 6.099 6.0995 6.1 6.1005 6.101 6.1015 6.102 6.1025 


35 


6.1085 6.109 6.1095 6.11 6.1105 6.111 6.1115 6.112 6.1125 


33 


6.1085 6.109 6.1095 6.11 6.1105 6.111 6.1115 6.112 6.1125 


30 


6.1285 6.129 6.1295 6.13 6.1305 6.131 6.1315 6.132 6.1325 
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TABLE HID 



Temperature 


PH 

7.34 7.33 7.32 7,31 7.3 7.29 7.28 7.27 7.26 


40 


6.083 6.0835 6.084 6.0845 6.085 6.0855 6.086 6.0865 6.087 


39 


6.083 6.0835 6.084 6.0845 6.085 6.0855 6.086 6.0865 6.087 


38 


6.093 6.0935 6.094 6.0945 6.095 6.0955 6.096 6.0965 6.097 


37 


6.093 6.0935 6.094 6.0945 6.095 6.0955 6.096 6.0965 6.097 


36 


6.103 6.1035 6.104 6.1045 6.105 6.1055 6.106 6.1065 6.107 


35 


6.113 6.1135 6.114 6.1145 6.115 6.1155 6.116 6.1165 6.117 


33 


6.113 6.1135 6.114 6.1145 6.115 6.1155 6.116 6.1165 6.117 


30 


6.133 6.1335 6.134 6.1345 6.135 6.1355 6.136 6.1365 6.137 



TABLE HIE 



Temperature 


PH 

7.25 7.24 7.23 7.22 7.21 7.2 


40 


6.0875 6.088 6.0885 6.089 6.0895 6.09 


39 


6.0875 6.088 6.0885 6.089 6.0895 6.09 


38 


6.0975 6.098 6.0985 6.099 6.0995 6.1 


37 


6.0975 6.098 6.0985 6.099 6.0995 6.1 


36 


6.1075 6,108 6.1085 6.109 6.1095 6.11 


35 


6.1175 6.118 6.1185 6.119 6.1195 6.12 


33 


6.1175 6.118 6.1185 6.119 6.1195 6.12 


30 


6.1175 6.118 6.1185 6.119 6.1195 6.12 



Since [TCO2] is very close to the sum of [HCO3'] and APCO2, [HCO3"] = [TCO2] - 
A.PCO2, then: 

Eq. 96 pH = pK + log(([TC02] - XPCO2)/ XPCO2) 

The degree of shift of the HODC can be calculated by: 

Eq. 97 temperature factor = antilog{0.024*(37-temperature)} 

Eq. 98 pH factor = antilog{0.48*(pH-7.40)} 

Eq. 99 base excess factor = antilog{-0.0013*base excess} 

Calculation of blood oxygen content can be made as follows: 

Eq. 100 content (ml O2 /dl) = THb(g/dl)*502* 1 .38(ml 02/g HbO2)+0.003*PO2 

Conversion of PO2 to SO2 (and of SO2 to PO2) is done using modifications of 
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Adair's equation or Kelman's computation: 

Eq. 101 SO2 = (25*(0.0257*PO2+2*0.00075*PO2^+3*0.00000444*PO2^ 

+ 4*0.00000255*PO2V(l + 0.0257*PO2 +0.00078*PO2^ 
+ 0.00000444*PO2^+0.00000255*PO2^)) 
Calculation of base excess can be done by the following formula, or other known means: 

Eq. 102 base excess = (1 - 0.0143*Hb)*([HCO3 ] - 24) 

These calculations and conversions represent approximate values as understood by 
those of skill in the art and the invention is not limited to them. As improved estimations 
become available, they may be used with the methods of the invention. Accordingly, 
alterations of this algorithm will be suggested to those skilled in the art, and are meant to 
be included within the scope and spirit of this application. 

Alternatively, values for PCO2 and bicarbonate can be computed without the use of 
temperature perturbation if the pH is known at two or more values. For example, 
techniques described elsewhere in this application show that we can perform simultaneous 
measurements for arterial and venous blood. The two compartments have the same 
temperature but different pH values, and this allows formation of the Henderson- 
Hasselbach formulation into two equations for two unknowns: PCO2 and [HCOs']. 

(D-6). Determination of Temperature 

The shifts in the HODC shown in Figure 39 relative to the curve shown in Figure 
38 are caused primarily by two factors: pH and temperature. Thus, the pH of two samples 
of blood at the same temperature will define which curve the blood is at, and thus what the 
temperature of the blood is. The techniques described elsewhere in this application can be 
used to provide simultaneous determinations of arterial and venous blood pH within a 
probe. As these pH values will be slightly different, the HODC can be defined, and 
temperature calculated quite accurately. This method of the invention gives a true, 
noninvasive measurement of blood temperature. Since the blood flow to the finger or 
other member within the probe is continuously refreshed, the measurement is a very 
accurate estimate of central or "core" temperature. 

(D-7). Determination of Glucose Concentration 

Most substances in blood are present in essentially the same concentration in 
arterial and venous blood. However, this is not true for substances that are consumed in 
significant quantities in the metabolic process. The most prominent of these are oxygen 
and glucose. This is especially true when dealing with molar quantities, as in 
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spectrophotometry. The technique for determining venous oxygen saturation h£is been 
described above. Once both arterial and venous saturation are determined, their 
absorbances can be "normahzed". This can also be performed for pH, as above. A 
subtraction of hemoglobin absorbance may then be performed, effectively removing it. 
What remains is the absorbance due to non-hemoglobin absorbers. This provides an 
estimation of glucose concentration since the largest non-hemoglobin absorber is glucose. 
Additionally, the absorbance at each wavelength can be compared for arterial and venous 
blood. A vast majority of the non-hemoglobin absorbers, but not glucose, will be of the 
same concentration in both arterial and venous blood. By subtracting out the constant 
absorbance between arterial and venous blood, the remaining glucose absorbance can be 
measured, and glucose concentration can be accurately determined. 
(D-8). Determination of Chemical Analvtes 

Once the concentration of hemoglobin is known, the concentration of many 
substances in blood could be determined as ratios of the concentration of hemoglobin. 
The use of multiple wavelengths will allow differentiation between molecules, and 
comparison can be made to concentration of hemoglobin even without knowledge of path 
length. A further embodiment of the invention establishes standards by creating layers or 
films containing known amounts of analytes of interest. These can then be compared 
optically with the changes due to blood volume shifts within a probe. Such techniques are 
not possible with conventional pulse oximetry as the signal from the film would simply be 
filtered out with all other unchanging (DC) signal. This invention overcomes this 
drawback by using venous measurements, and thus does not require limitation to pulsatile 
(AC) signal. 

One exemplary method includes the steps of placing a film that contains a known 
concentration of an analyte within a probe system, such as one configured for a finger. 
The absorbance of the finger plus film system is measured for a given number of 
wavelengths of light. The blood volume of the finger is changed by raising or lowering it 
with respect to the heart (as outlined elsewhere in this application). As it can be assumed 
that the change in absorb£mce is due to change in blood volume, one can compare blood 
absorbance to film absorbance, and hence relative concentration of the analyte in blood to 
that in the film. Sufficient wavelengths must be examined to provide specificity for the 
analyte in question, as there are many potential absorbers in blood. 
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(D-9). Determination of Congenital Heart Disease and Anatomic Anomalies 
Diagnosis of many disorders with anatomic anomalies can be made by the 
detection of unexpected propagation times, and abnormal propagation delays between 
right- and left-sided organs. The ability to measure both arterial and venous saturation, as 
5 well as arterial and venous pressures, can aid ftirther in investigations. 

For example, coarctation of the aorta is a common congenital abnormality. It is a 
narrowing or constriction of the thoracic aorta. It can occur before or after the origin of 
the left subclavian artery. In either case, the blood pressure will be much lower in the legs 
than the arms. Because of the constriction, both pressure and flow waves will be greatly 
10 delayed in the lower limbs as compared to the upper. If a child has a coarctation proximal 
to the origin of the left subclavian artery, there will be a significant discrepancy in the 
pulse and flow wave arrival times of the two arms, with the right preceding the left. 
, ^ (D- 1 0) . Determination of Dysrhythmias 

□ By measuring blood pressure and the electrocardiogram simultaneously, the 

P5 15 diagnosis of dysrhythmias can be aided greatly. Both arterial and venous pressure are 
recorded simultaneously with the ECG, allowing differentiation of atrial vs. ventricular 
arrhythmias as well as correlation to Q-T interval. 

It can be difficult to distinguish dysrhythmias of atrial conduction from those of 
r- ventricular conduction. However, these often manifest much different hemodynamic 

1^220 effects. Those with normal ventricular conduction and eject often maintain cardiac output 

and blood pressure closer to normal. 
M (D- 1 1 V Determination of Additional Cardiovascular Parameters 

By measuring blood pressure and the electrocardiogram simultaneously, many 
additional parameters, such as systolic and diastolic pressure time indices, can be 
25 determined. An enormous amount of information can be gleaned from the use of probes 

on opposite sides of the body combined with hydrostatic perturbations. The time of arrival 
of a pulse to paired members is different, but the velocity of the pulse is also different. 
Examination of pulse propagation time, pulse propagation phase or delay, pulse velocity, 
and pulse amplitude gives four parameters that may change in different ways for each 
30 perturbation. Circuitry such as that shown in Figure 13, is suitable for use in such 

embodiments. In one embodiment, raising and lowering an arm by the same amount may 
give different changes. Raising and lowering the other arm by the same amount may give 
still different changes. In an alternative embodiment, raising an arm by a given amount. 
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then raising again by the same amount, may give different changes. Raising the other arm 
by the given amount, then raising again by the same amount, may give still different 
changes. Similarly measurement can be taken for lowering. 
(D- 1 2\ Determination of Aortic Pressure 

Using the methods and probes of the invention allows the aortic pressure to be 
determine and the reconstruction of the aortic pressure curve. This can be accomplished 

through fluid mechanics or changing the reflection wave from the lower extremities. 
Details of exemplary techniques may be found in United States Patent Nos. 5,265,011, 
5,882,3 1 1 , and 6,010,457 which are hereby incorporated by reference thereto. 
(D-13y Determination of Carotid Stenosis 

The carotids are paired euteries running from the chest through the neck to the 
head, and provide the major source of blood supply to the head and brain. As they have 
different origins, they have different path lengths for pulse transmission and blood flow 
from the heart. A significant discrepancy in pulse and flow wave velocity in one carotid 
artery as compared to the other, or in paired arteries arising from the carotids (such as the 
facial arteries) suggest a stenosis on the side exhibiting the delay. 

(D- 1 4). Determination of Peripheral Vascular Disease 

Peripheral vascular disease results in narrowing of the lumen of peripheral arteries, 
most commonly the large arteries of the leg. Because of this stenosis, the blood flow 
characteristics of the artery affected are altered, and thus the pulse and flow waves distal 
to the obstruction. If a probe is placed distal to the obstruction, significant delays will be 
noted. Again, this is best done when compared to a paired organ; for example, placing 
probes on the great toes of both feet to compare the wave characteristics of the arteries of 
the two legs. 

(D-15y Determination of Compartment Syndrome 

Compartment syndrome is hypoperfusion and resulting ischemia of a body 
member, commonly an arm or a leg, secondary to swelling caused by tissue damage. If 
the swelling is significant, the increased volume of the tissue may increase the pressure in 
the limb (compartment) to such an extent that it cuts off the arterial blood supply. This 
can be monitored by placing a probe distal to the compzirtment to measure the pulse and 
flow wave velocities. This is done preferably with a paired organ to compare the two. 
Slowing of the pulse wave and increase in the arrival time signal diminishing blood 
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pressure within the compartment, which can eventually resuh in ischemia and gangrene if 

steps are not taken to relieve the prsesure, normally with surgery. 

(D- 1 6). Determination of Autonomic Activity, Sedation and Depth of Anesthesia 
There are normal vziriations in pulse baseline, pulse amplitude, and beat-to-beat 

variability. These factors are thought to be indicative of autonomic nervous activity or 

"tone". All these factors can be measured by the device and techniques of the invention. 

An alternative embodiment invention cein also be used as a monitor of autonomic activity, 

sedation, or "depth of anesthesia". 

(D- 1 7). Determination of Heart Rate 

According to the invention, heart rate can be determined by counting the pulsatile 
arterial signal for a known length of time, or by the ECG impulse. 
(D- 1 8\ Determination of Respiratory Rate 

The impedance changes of the chest due to filling and emptying can be measured 
from the electrocardiogram tracing. During normal breathing, negative pressure is created 
within the chest by lowering of the diaphragm and expansion of the rib cage. This 
negative pressure causes blood to empty more rapidly from the peripheral into the central 
veins. This is also the case when respiration is assisted by a negative-pressure device such 
as the "iron lung". 

During modem mechanically-assisted ventilation (with "ventilators"), positive 
pressure is created within the chest by forcing air into the lungs. For both positive- and 
negative-pressure ventilation, expiration is passive. This respiratory veiriation by itself can 
be used as an estimate of cardiac Ailing, giving left heart pressures. This determination 
can be assisted by the use of the hydrostatic techniques described above. 

fD-191. Determination of EEG and EMG 

Just as leads capable of detecting and transmitting electrocardiogram (ECG) 
signals can be included in the probes (as outlines elsewhere in this application), similar 
techniques will allow monitoring of the electroencephalogram (EEG) and electromyogram 
(EMG). 

(D-20). Determination of Cardiac Ejection Fraction 

This determination is made based upon the rate of change of saturation. If blood in 
a container is of a stable, known saturation, and blood of a different saturation is added, 
the rate of change is dependent on two things: the saturation of the blood being added. 
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and the amount of blood that is flowing in (equivalent to blood flowing out). The stroke 
volumes of the right and left ventricles must be the same, so the inflow per heartbeat to the 
left ventricle must be the same as the outflow. While the saturation of blood flowing from 
the lungs to the left heart may not be precisely known, it will quickly equilibrate. This 
5 equilibration process allows the determination of ejection fraction to be made. 

The saturation of the blood in the left side of the heart, specifically the left 
ventricle, is completely dependent on the oxygenation of the blood pumped into the lung 
by the right ventricle. The blood pumped by the left ventricle into the arterial tree 
maintains its saturation until it reaches the capillary beds of the organs. The blood 
10 returning to the right heart is a combination of widely varying saturations (due to the 

greatly differing extractions of the various organs). However, the blood pumped by the 
right ventricle equilibrates over the vast capillary network of the lungs. Furthermore, any 
differences in saturation of blood in the pulmonary venous system are thoroughly mixed in 
p the left atrium and left ventricle. 

p 15 Breathing an oxygen-enriched gas mixture will change the alveolar concentration 

M of oxygen. A single deep breath after exhalation may be the preferred method. While 

healthy subjects can inhale a vital capacity (VCII) breath of near total lung capacity 

in 

^ * (TLC), one cannot expect this of all patients. However, it is not necessary to know the 

M precise oxygen level in the lungs. All that is necessary is to induce a change in saturation, 

f 2 20 This works well for the determination, as it is much easier to produce saturation change in 
^ disease states, because the patient starts farther away from the maximum value of 100%. 

M Preferably, this may be achieved through a breath-hold after the inspiration. Wash-out can 

be accounted for (one knows that the inspired gas is air), however, this complicates the 
otherwise simple mathematics. The breath holding is required for only several heartbeats, 
25 normally only a few seconds. 

Even with holding one's breath, oxygen levels within the lung drop due to 
absorption by blood. However, even with high hemoglobin levels, the normal amount of 
oxygen carried by blood is about 20 cc/dL, and, as mixed venous saturation is normally 
about 75%, only about 5 cc/dL are taken up. Given a cardiac stroke volume of 1 00 cc 
30 (higher than normal), this is only 5 cc per heartbeat. A small breath will provide enough 
oxygen for several heartbeats, which should be all that is required. 

Thus, an approximation is made that the concentration of oxygen in the lungs after 
oxygen inhalation is stable for several heartbeats. The first heartbeat after the inhalation 
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will result in blood of a changed saturation entering the left heart from the lungs (maybe 
2nd heartbeat, accounting for "dead space" of the pulmonary veins). This will mix with 
the blood remaining from the last cardiac ejection. This remaining blood is the end- 
systolic volume of the left ventricle. It is the end-diastolic volume minus the stroke 
volume. The saturation of the blood in the left ventricle after mixing is between that of the 
end-systolic blood and the fresh blood from the lungs. The saturation of the mixed blood 
is measured after it is ejected to the periphery. The arrival of the blood at the periphery 
can be judged by the change in saturation from the baseline measurement. The saturation 
rises depending on the ratio of end-systolic blood to stroke volume. The higher the stroke 
volume, the faster the change in saturation. This can be looked at as either a "wash in" or 
"wash out" function. 

Thus, a high ejection fraction results in a rapid change in saturation as shown in 
Figure 40. The rapid change is represented by the relatively steep slope of the plot. In 
contrast, a low ejection fraction causes a slow change in saturation as shown in Figure 41, 
Correspondingly, the slow change is represented by a relatively shallow slope of the plot. 
The exact same physiologic process reveals information about the cardiac output. As 
blood leaves the left ventricle, it travels to the aorta and then to the periphery. Just as the 
pulse arrives at different points of the body at different times, so does the blood flow from 
a given heartbeat. The blood travels much more slowly than the pulse wave, so the time 
delays are much greater. The change in oxygen saturation produces identifies the blood 
from individual beats to trace its arrival time to different parts of the body. By placing 
pulse oximeter probes at different points of the body (a finger and a toe, for example) and 
continuously monitoring saturation, the saturation change produced by the oxygen 
inhalation will occur in all parts of the arterial tree, but at different time intervals. The 
timing and ratio of these intervals is determined by the cardiac output - the longer the 
arrival time of the saturation change (and the greater the interval between finger and toe), 
the lower the cardiac output. Conversely, the shorter the arrival time of the saturation 
change (and the shorter the interval between finger and toe), the higher the cardiac output. 
Note that it is not necessary to maintain a breath-hold through the measurement since only 
a change in saturation is required. 

An exemplary determination of cardiac ejection fraction starts with baseline 
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saturation and saturation increasing because of inhalation of oxygen-rich gas: 

• Let dS be change in saturation S. 

• After 1 time constant, dS will have risen to approximately 63% of its final value. 

• After 2 time constants, dS will have risen to approximately 86.5% of its final value, 

• After 3 time constants, dS will have risen to approximately 95% of its final value. 

• After 5 time constants, dS will have risen to approximately 99% of its final value. 

• If the resting ("control") saturation is 95%, and it rises to 99% after oxygenation, total dS 
is 4.0%. Thus, 1 time constant is 0.63 * 4.0 = 2.52 percent, for a saturation value of 95.0 + 
2.5 = 97.5%. 

• 2 time constants are 0.865 * 4,0 = 3.46 percent, for a saturation value of 95.0 + 3.46 
= 98.5%. 

• 3 time constants are 0.95 * 4.0 = 3.80 percent, for a saturation value of 95.0 + 3. 
= 98.8%. 

• 5 time constants are 0.99 * 4.0 = 3.96 percent, for a saturation value of 95.0 + 2. 
= 98.96%. 

For this system, the time constant is the number of heartbeats that it would take to 
fill the ventricle from the empty state, or empty the ventricle from the full state. For 
example, an ejection fraction of 50% implies two heartbeats to fill or empty the ventricle, 
so the time constant is 2 beats. An EF of 80% implies a time constant of 1 .25 beats, 
although this might be expressed more physiologically as 5 heartbeats are 4 time 
constzints. A "worst case" scenario of an EF around 20% corresponds to a time constant of 
5 heartbeats. Clinically, some type of calibration may be required, as actual measurements 
may not correspond exactly to these theoretical values. 

It is expected that a minimum of 2 time constants are required to produce a 
satisfactory measurement, so an EF 20% would take 10 or more heartbeats. However, 
there is the trade-off in that there are more beats to analyze, so the curve can be fitted more 
precisely. Thus, it should be possible in almost all cases to produce an accurate 
determination within a few heartbeats after noting the saturation change. 

Without departing from the spirit and scope of this invention, one' of ordinary skill 

can make various changes and modifications to the invention to adapt it to various usages 
and conditions. As such, these changes and modifications are properly, equitably, and 
intended to be, within the full range of equivalence of the following claims. 
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